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Abstract 

The paper contains an exposition of recent as well as sufficiently 
old results on determinantal random point fields . We start with 
some general theorems including the proofs of the necessary and suf- 
ficient condition for the existence of determinantal random point field 
with Hermitian kernel and a criterion for the weak convergence of its 
distribution. In the second section we proceed with the examples of 
the determinantal random fields from Quantum Mechanics, Statistical 
Mechanics, Random Matrix Theory, Probability Theory, Representa- 
tion Theory and Ergodic Theory. In connection with the Theory of 
Renewal Processes we characterize all Hermitian determinantal ran- 
dom point fields in R 1 and Z 1 with independent identically distributed 
spacings. In the third section we study the translation invariant de- 
terminantal random point fields and prove the mixing property of any 
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multiplicity and the absolute continuity of the spectra. In the last sec- 
tion we discuss the proofs of the the Central Limit Theorem for the 
number of particles in the growing box and the Functional Central 
Limit Theorem for the empirical distribution function of spacings. 

1 Definition and General Properties of De- 
terminantal Random Point Fields 

Let E be a one-particle space and X a space of finite or countable configu- 
rations of particles in E. In general E can be a separable Hausdorff space, 
however for our purposes it is enough to consider 

m 

E = Y[Ej, where Ej = R d (or Z d ) (1.1) 

j'=i 

If it is not mentioned specifically otherwise we always assume below E = M. d 
with the understanding that all resuls can be easily generalized to (1.1). We 
assume that each configuration £ = (xi),Xi G E,i G Z 1 (or Z\ if d > 1), is 
locally finite, that is for every compact K C E = G K) is finite. 

The particles in £ are ordered in some natural way, e.g., X{ < for d—1, 
and if d > 1 then either Xi = x i+i , or 

where Xi = {xf \ . . . ,x^), or \xi\ = \xi + \\ and there exists 1 < r < d such 

(i) (i) (r) (r) 

that x[ < xY+i, 1 < J < r — 1, and < 

To define a a-algebra of measurable subsets of X we first construct the 
so-called cylinder sets. Let B C E be any bounded Borel set and n > 0. We 
call C,f = {£ G X : #b(£) = n} a cylinder set. We define £> as a cr-algebra 
generated by all cylinder sets (i.e., B is a minimal a-algebra that contains all 

Definition 1. A random point field is a triplet (X,B,P) where P is a 
probability measure on (X, B) . 

This definition raises a natural question, namely how one can construct 
such probability measures. The corresponding theory was developed by 
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Lenard in [L1-L3] where a general case of E locally compact Hausdorff space 
satisying the second axiom of countability was studied. If E = IR d or 7L d one 
can proceed quite naively by employing Kolmogorov's fundamental theorem 
from the theory of stochastic processes ([K]). Let t and s be two vectors from 
E with rational coordinates t = (t^, . . . , t^), s = (s^, . . . , s^). We denote 
an open rectangle {x = (x w , . G E : x ij) = + 6j(s® - t (j) ),0 < 

0j < l,j = 1, . . .,d} by \l tjS . Let us denote the family of finite unions of 
open, closed or semi-closed rectangles with rational t, s by 1Z. Suppose we 
are able to construct a joint distribution of non-negative integer- valued ran- 
dom variables t}d,D G 1Z (that we later identify with such that the 
following finite-additivity condition holds 

n 

Vd = ^riDt (a.e.) (1.3) 

i=i 

if D = |J" =1 Di, D, Di e TZ, % — 1, . . . n. One immediately can replace (1.3) 
then by (j-additivity property 

oo 

Vd = ^VDi (a.e.), (1.4) 

i=i 

oo 

D = I I Di, D,Di G 1Z, i = 1,..., (of course the fact that i] D takes only 
i=i 

non-negative integers is essential here!). 

It is then easy to see that the joint distribution of random variables 
#d = Vd, D G TZ with (1.3) (or (1.4) for that matter) uniquely defines a 
probability distribution on (X,B). 

Since in many cases it is convenient to define the distribution of random 
variables through their moments the following definition appears natural: 

Definition 2. Locally integrable function pk : E k — > is called the /c-point 
correlation function of the random point field (X, £>, P) if for any disjoint 

TO 

bounded Borel subset A ± , . . . , A m of E and hi G Z\, i — 1, . . . m, ^ h L = k 

i=i 

the following identity holds: 

E 7H VTT = / p k {x 1 ,...,x k )dx 1 ...dx k (1.5) 
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where by E we denote the mathematical expectation with respect to P. In 
particular pi(x) is the density of particles, since 



E#a = / pi(x)dx 



A 



for any bounded Borel A C E. In general p k (x 1: . . . ,x k ) has the following 
probabilistic interpretation: let [x±,Xi + dxi], i = 1, . . . , k be infinitesimally 
small boxes around Xi, then Pk(xi,X2, ■ ■ ■ ,x k )dxi ■ . . . ■ dxk is the probabil- 
ity to find a particle in each of these boxes. The problem of existence and 
uniqueness of a random point field defined by its correlation functions was 
studied in [L1-L3]. Not very surprisingly, Lenard's papers revealed many 
similarities to the classical moment problem ([A], [S2]). In particular the 
random point field is uniquely defined by its correlation functions if the dis- 
tribution of random variables is uniquely determined by its moments. 
The sufficient condition for the uniqueness derived in [LI] reads 



E 

fc=0 



(jj—r x, / Pk+j(xi,...,x k+j )dx 1 ,...dx k+j ) =00 (1.6) 



for any bounded Borel A C E and any integer j > 0, however we invite the 
reader to check that the divergence of the series with j = 0, namely 



00 

E 

k=0 



J ^Pk(xi,...,x k )dx 1 ,...dx k ^j =00 (1.6') 



implies (1.6) for any j > 0. In [L2], [L3] Lenard obtained the necessary 
and sufficient condition for the existence of a random point field with the 
prescribed correlation functions. 

Theorem 1. (Lenard) 

Locally integrable functions p k : E k — ► M, 1 , k = 1,2, . . . are the correla- 
tion functions of some random point field if and only if the Symmetry and 
Positivity Conditions below are satisfied. 

a) Symmetry Condition 

p k is invariant under the action of the symmetric group S k , i.e., 

p fe (aV(i), • • • ,x a{k) ) = p k (xi, ...,x k ) (1.7) 

for any a G S k . 
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b) Positivity Condition 

For any finite set of measurable bounded functions (f k : E k — > M 1 , 
k — 0, 1, . . . , N with compact support, such that 

N 

¥o + Y^ Y ( Pk(x il ,...,x ik )>0 (1.8) 

fc=l i 1 ^---y^i k 

for all £ = (xi) G X, the next inequality must be valid: 

N 

Vo + y^ Vk(xi, x k )p k (xi, x k )dxi ...dx k >0. (1.9) 
k=i JeU 

The necessary part of the theorem is quite easy since both conditions have 
an obvious probabilistic interpretation. In particular the Positivity Condition 
means that the mathematical expectations of a certain class of non-negative 
random variables must be non-negative. The sufficient part is more elaborate 
and relies on an analogue of the Riesz Representation Theorem and the Riesz- 
Krein Extension Theorem (a close relative of the Hahn-Banach Theorem). It 
should be noted that Lenard established his results in a general setting when 
E is locally compact Hausdorff space with the second axiom of countability. 

One can obtain a slightly weaker (but still hopelessly ineffective!) variant 
of the Positivity Condition by approximating ip k from above by step func- 
tions. Let V k be the class of polynomials in k variables that take non-negative 

k mi — 1 

values on non-negative integers. Since the polynomials {J^J | J j), mi > 

i=l j=Q 

0} form a linear basis in the vector space of all polynomials in k variables, 
we can represent any q(x±, . . . , x k ) G V k as 

k rrii — l 

q(x 1 ,...X k )= Yl a mi,..,m k ■ IJ IJ (Xi - j) (1.10) 

m.i,...,mfc>0 i=l j=0 



Positivity Condition*: For any q G Vk,k > 1, any bounded Borel sets 
Ai, . . . , A k C E, the following condition must be satisfied: 

°0,...,0 + / , / J a m\,...,m k I Pm(Xi, . . . , X m )dXi . . . dx m > 0. 

— — JT\ k A mi 

m>lmH hmj.=m ili=l^i 

(1.11) 
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Indeed, the l.h.s. at (1.11) is equal to 

Eq(# Al ,...,# Ak ) =E[a ,... + ^ Yl a ™u-,™k ' Yl 

m>imH Ym k =m i^—^i m (1.12) 

One can notice that in a sense the Positivity Condition* is similar to the 
condition on the moments of the integer-valued nonnegative random variable. 

In our paper we will study a special class of random point fields introduced 
by Macchi in [Ma] (see also [DVJ]). We start with an integral operator K : 
L 2 {K d ) — > L 2 (M. d ) that we assume to be non- negative and locally trace class. 
The last condition means that for any compact 5cl d the operator K ■ xb 
is trace class, where Xb(x) is an indicator of B. Therefore we have 

K>0, Tt( X bK ■ xb) < +00 (1.13) 

The kernel of K is defined up to a set of measure zero in ~K d x W 1 . For our 
purposes it is convenient to choose it in such a way that for any bounded 
measurable B and any positive integer n 

Tt(( Xb Kxb))= [ K(x,x)dx (1.14) 
Jb 

It appears that one can indeed achieve this. We start with 

Lemma 1. ([S3], [AvSS] Remark 3.4) Let K be trace class on L 2 (R d ). Then 
its integral kernel may be chosen so that the function M(x, y) = K(x, x + y) 
is a continuous function of y with values in L 1 (M. d ). Furthermore if m(y) = 
J M(x,y)dx, then TrK = m(0) = J K(x,x)dx. 

Proof. We give the proof only when K is non-negative. The general case is 
quite similar. Let {\j}j>\ is the set of non-zero eigenvalues of K and {<Pj}j>i 
is the set of the corresponding eigenfunctions. The canonical form of K (as 
a self adjoint compact operator) is 

K = Y, x r(<Pi>')'<Pi (1-15) 



Fix y G M d and consider M(x, y) = Xj ■ ifj(x) ■ ipj(x + y) as a function of 

3=1 

x. Since ||^-(-) • ipj(- + y)\\i = j Rd \(pj(x) ■ ipj(x + y)\dx < ||^|| 2 • ||^|| 2 = 1, 
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the series defining M(-,y) converges in L 1 (IR d ) for any y and ||M(-,y)|| 1 < 

oo 

we now consider K(x,y) = M(x,y — x), it is well 

3=1 

defined for a.e. (x, y) G l rf xR d and gives a kernel for K. The /^-continuity of 

N 

M (•, y) follows from || ^ A, (yj(-)<fj(- + yi) - <Pj(-) ■ <fj(- + y 2 )) II 1 < YI V 

3>1 3=1 

W^PjlU ' W^Pji' + Vi) ~ + l/2)l|2 + Aj. Choosing TV sufficiently large one 
can make Aj < |. Choosing y x sufficiently close to y 2 so that for each 

j>N 

1 < j < ^ — + 1/2 — yi) || 2 < tttI — r- we have the first term less 

2'2^3 = 1 A) 

than I as well. □ 
With the help of Lemma 1 we derive 

Lemma 2. Let K be a non-negative locally trace class operator on L 2 (M d ). 
Then its integral kernel can be chosen in such a way that for any bounded 
measurable B C ~R d the function 

M B (x,y) = (K ■ xb)(x,x + y) 

is a continuous function of y with values in L l (B). Furthermore, 

Tt(XbK X b) = [ K(x,x)dx 

JB 

Proof. Let K n = X[-n,n] d ' K ' X[-n,n] d - By Lemma 1 one can choose the kernel 
K n (x, y), such that K n (-, -+y) is continuous in n, n] d )-norm. We denote 
M n (x,y) = K n (x,x + y). Since K n+1 (x,y) = K n (x,y) for almost all (x,y) G 
[— n, n] d x [— n, n] d we conclude that for almost all |y| < n M n+1 (x,y) = 
M n (x,y) for a.e. \x\ < n— \y\. The /^-continuity of M n+ i(-, y), M n (-, y) 
allows to replace "almost all \y\ < n" by "all \y\ < n" . Therefore for any y 
the values of M n (x,y) eventually agree for a.e. x. We denote this value by 
M(x,y). The function M(-,y) inherits local /^-continuity from {M n (-,y)} 
and (1.14) follows. It may be worthwhile to note that for any positive integer 
k and bounded measurable B 1 C M d , ... , B k C W 1 the kernel 

M% ] (x, y)=((K- X B 1 )-...-(K- X B k )) (x, x + y) (1.16) 

< — fctimcs > 
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is also a continuous function of y with the values in L 1 (R d ) and 

Tr(K ■ xbJ • • • • • {K ■ XB k ) = / K(x 1 ,x 2 ) ■ ■■■■ K(x k , x^dxi . . . dx k , 

J Bix — xB k 

(1.17) 

in particular 

TriK ■ xb) ■ ■ ■ ■ ■ {K ■ xb) = / K(x 1 ,x 2 ) ■ ■■■■ K(x k , x 1 )dx l . . . dx k . 

JB k 

(1.18) 

Indeed, for any two versions of the integral kernel and k > 1 the expressions 
in the last two integrals coincide up to a set of measure zero. Since we have 
already proved that there exists the kernel of 

K-XB x ...K- X B k (1.19) 

satisfying the L 1 - continuity condition from Lemma 2, the same condition is 
satisfied by any variant of the kernel. 

□ 

Definition 3. A random point field in E is called determinantal (or fermion) 
if its n-point correlation functions are given by 

p n (x 1 ,...,x n ) =det(iir(a:i,a: J -)) 1 < i < n (1-20) 

In the case E = \jf =1 Ej,Ej = R d the definit ion takes the following form: 
Let K be a trace class operator on L 2 (R d ) © • • • © L 2 (R d ). Then K has a 

< m times > 

matrix valued kernel (K rs (x, y))i< r , s <m, x, y G R d . 

Definition 3'. A random point field in E is called determinantal (or fermion) 
if its n-point correlation functions are given by 

Pni^lli 2-12? • • • ; • • • ; %m2i ■ ■ ■ i %mim) ^et {K rs (x r j , 3^aj))l<i<i I . J r=l,...,m) 

l<j<i s ,s— 1, m 

(1.21) 

where n = %\ + i 2 + ■ ■ ■ + i m , x ri G E r , 1 < r < m, 1 < i < i r . 

Remark 1. If the kernel is Hermitian-symmetric then the non- negativity of 
n-point correlation functions implies that the kernel K(x,y) is non-negative 
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definite and therefore indeed K must be a non-negative operator. It should 
be noted however that there exist determinantal random point fields cor- 
responding to non-Hermitian kernels (see the remark after (1.36) and the 
examples in the sections 2.2 and 2.5). 

Remark 2. The condition (1.13) is satisfied for all continuous non-negative 
definite kernels (see [GK], section III. 10 or [RS], vol. Ill, section XI.4). In 
general situation when K(x, x) is locally integrable, non-negative definiteness 
of K(x, y) implies that Kb is a Hilbert-Schmidt operator and one can use a 
theorem of Gohberg-Krein ([GK], section III. 10, theorem 10.1) that claims 
that a non-negative Hilbert-Schmidt operator A is trace class iff 

I r d 

limpp J J\j2h - \x> - yJ\] + A(x,y)dxdy < oo (1.22) 

where t+ = max(t,0), x = (x 1 , . . . ,x d ), y = (y 1 ,...y d ), and TrA is then 
given by (1.22). 

An interesting generalization of determinantal random point fields, so 
called immanantal random point fields (processes) was introduced by Diaco- 
nis and Evans in [DE]. 

The classical formula of Fredholm (see [SI], Chapter 3) claims that a 
trace class operator with a continuous (in a usual sense) kernel satisfies 

Tr(A"(A)) = iy det(A(x i ,x j )) 1 ^ n dx 1: ...dx n (1.23) 

In general the kernel K(x,y) may not be continuous, however (1.18) and the 
Lidskii theorem (see e.g., [RS], volume IV, section XIII. 17 or [SI], Theorem 
3.7) imply 



/ 



K(x\, x 2 ) ■ ... ■ K(x n , Xl )d Xl ...dx n = Y, WKb), (1-24) 

3=1 



Tr(A n (K B ))= J2 K{Kb)--..-K(Kb) (1-25) 

31<-<3n 

Combining (1.24) and (1-25) one arrives at 

Tr(^{K B )) = - f det^,^))^ dx u ...dx n (1.26) 

n\ J fin - ' J - 
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It follows then from (1.17) that 



Tr((K.x Bl )A---A(K. X B n )) = 

1 f (1-27) 
— / det (K(xi, Xj) ■ Xb 3 (xj)) 1 < ij < n ^i ■■■dx n 

Definition 4. Let the kernel K as in Lemma 2. We say that it defines a 
determinantal random point field (X,B,P) if (1.21) holds. 

Theorem 2. Let (X, B,P) be a determinantal random point field with 
the kernel K . For any finite number of disjoint bounded Borel sets Bj C 
E, j = l,...,n, the generating function of the probability distribution of 
#b j = e Bj} is given by 

n n 

E U z f Bj " =det(Id + XB^2(z j -l)-K- XBj ) (1-28) 
j=i j"=i 

Remark 3. (1.28) is the equality of two entire functions. The r.h.s. of (1.28) 
is well defined as a Fredholm determinant of a trace class operator (see e.g., 
[RS], volume IV, section XIII. 17 or [SI], section 3). 
Recall that by definition 

n oo n 

E H z#Bj= E P(#Bi = k j ,j = l,...,n)-l[zj i (1-29) 

i=i fci,...,fe„=o j=i 

and 

oo n m 

(1.30) 



detfld + XB^-^-^-XB, j =1 + E E II 

^ j=l ' m=lj u ...j m =ie=l 

fat - !) • Tr (Xs -K-xb^A-'-A XbK ■ X B jm ) 



Proof of Theorem 2. The Taylor expansion of the generating function 
near (z u . . . , z n ) = (1, . . . , 1) is given by 

n oo n / ii \ i n 

eTT^=i + E E eTT ^ m -II^- 1 )^ 

11 J ^ ^ 11 (# B . - m,)! • (m,)! 11 VJ 7 

j=l m=lmH hm„=m j=l v " D J v Jy j=l 

(1.31) 
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The radius of convergence of (1.30) is infinite since 
Tr(iT • X B h A---AK- XBjm < —Tr(K ■ XB ) m , where B = [] B y (1.32) 

Therefore, it is enough to show that the coefficients in the series (1.30), (1.31) 
coincide. The case n = 1 follows then from (1.5), (1.21), (1.26). Using (1.27) 
instead of (1.26) we prove the case n > 1 as well. 

Remark 4. Theorem 2 is well known in the Theory of Random Point Fields 
(see [DVJ], p. 140, exercise 5.4.9) and in the Random Matrix Theory (see 
[TW1]). 

As we already mentioned above, if an operator K defines a determinantal 
random point field it must be non-negative because of the non-negativity of 
the correlation functions. It follows from Theorem 2, formula (1.28) that 
K must also be bounded from above by the identity operator, i.e., K < 1. 
Indeed, suppose > 1. Then there exists a bounded Borel B C E such 
that \\K B \\ > 1 + &i > i. Let Xi(K B ) > X 2 (K B ) > X 3 (K B ) > ... be the 
eigenvalues of K B and choose < zq < 1 so that 1 + (zq— 1)-Xi(K b ) = 0. Then 
Ezf B = Er=i p (#s = k ) z o = det(ld + (z -l)-K B ) = (by Theorem XIII. 106 
from [RS])= flj>i( l + ( z o - 1) • \{K B )) = 0. Therefore P(# B = k) = for 
any k, a contradiction. On the other side assume < K < 1 and let (1.28) 
define what we hope to be the distribution of non-negative integer-valued 
random variables 

Lemma 3. Let < K < 1 and K be a locally trace class operator. Then 
(1.28) defines the distribution of non-negative integer-valued random vari- 
ables with the additional property that for B = UiLi ^ 

n 

#B = ]T#B t (o.e.). (1-33) 

i=i 

□ 

We need to show three things: first, that (1.28) defines some finite- 
dimensional distributions; second, that the finite-dimensional distributions 
satisfy the additivity property (1.33); and third, that the finite-dimensional 
distributions are consistent and therefore we can apply the Kolmogorov's 
Fundamental Theorem to prove the existence of the distribution of 
Since the Fredholm determinant in (1.28) is 1 when Zi — 1, % — 1, . . . , n, 



11 



the first statement would follow from the non- negativity of the Taylor co- 
efficients of the Fredholm determinant at = 0,i = l,...,n. Consider 
< Zi < l,i = l,...,n and assume for a moment \\K\\ < 1 (the case 
Hi^ll = 1 would be treated later by a limiting argument). Let B = \_\™ =1 Bi. 
Then \\K B \\ < 1 and (Id — Kb)- 1 is a bounded linear operator such that 
(Id — Kb)- 1 — Id = K B ■ (Id — Kb)' 1 is trace class. Applying Theorem 
XIII.pl05 from [RS], vol. IV we obtain 



det(Id + X b - 1) • K ■ XBj ) = det((Id - K B ) ■ (Id + 

3=1 3=1 



Z 3 



n 



(Id - Kb)- 1 -Xb-K- xb 3 )) = det(Id - K B ) ■ det(Id + • 

(1.34) 



°3 

3=1 

oo n 

(Id - Kb)- 1 ■ XbK ■ xb,) = det(Id -K B )-J2 Tr (A fc (]T z, 

k=l 3=1 



(u-k b) -i. xb k. xb ,) )= y. (t \yif" )! n 

fcl,...,fc„>0 n ' 3=1 

z k / • det(Id - K B ) ■ TV(A" (A^( Xs ,(Id - T^)" 1 • • Xs .))) 



One can see from (1.34) that Taylor coefficients are, up to some positive 
factors, the traces of the exterior products of the non-negative operators, 
and, therefore, non-negative. We conclude that (1.28) defines some finite- 
dimensional distributions. Since E n™=i z#Bi = det(Id+x_B YH=i( z ~ ^)Kxb % ) 
det(Id+ (z - 1) • K B ) = Ez# B , we conclude that # B = f% =1 # Bl (a.e.). The 
formula (1.28) defined the finite dimensional distributions of # Bi for disjoint 
compact sets. In the case of non-empty self-intersections one represents Bi as 
UC^, where {Ck} are disjoint sets, defines distributions of #c k and then uses 
the additivity property (1.33) to define the distributions of To prove 

the consistency of the finite-dimensional distributions we note that (1.33) 
allows us to check it only for the disjoint B 1 , . . . , B n+ i. But then it trivially 
follows from det(Id + xb Y%=i( z j ~ X ) ' K ' Xb, ■+ Xb(1 - 1) • K ■ Xs n+1 ) = 
det(Id + XbY^=i( z j — 1) • K • Xb,)- The case \\K\\ < 1 is proven. Now let 

\\K\\ = 1. Denote by K^ := K ■ (1 — e),e>0 and the random variables 
corresponding to the kernel K^ e \ Since H-f^H < 1 the arguments above 
establish the result of Lemma 3 for K^ e \ It is an easy exercise to see that 
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^Yli=l Z i ' = 1 +E ra =l J2mi+-+m n =m^'Ilj=l (# (e) -mj)\-{mj)\ ' 11.7=1 ( Z l ~ V 



uniformly converges with all derivatives to E HILi ^f^* 011 compact sets as 

Lemma 3 is proven. □ 

The results above prove 

Theorem 3. Hermitian locally trace class operator K on L 2 (E) defines a de- 
terminantal random point field if and only ifO<K<l. If the corresponding 
random point field exists it is unique. 

□ The necessary and sufficient condition for the existence of the field has been 
already established. The uniqueness result easily follows from the general cri- 
terion (1.6') since ^ f Ak p k (xi,..., x k )dx 1 ...dx k = Ti(A k (K A )) < Tr( f, A) < 

i ' □ 

Consider arbitrary bounded Borel set B C E. Then Tr (Kb) = < 00 
and the number of particles in B is finite with probability 1. Let us write 
x = Uo<fe<oo C k' where as before Cf = {£ e X : #s(£) = k}. We choose a 
kernel for xb-K-xb in such a way (see Lemma 1) that (xb-K -xb)(x, x+y) = 
YllLi^B)^^) • (fii(x + y) is a continuous function of y in L l (B) norm. 
Assume for a moment that K B < 1. Then 

L B (x, x + y) = I _\.(b) ' ' Vi(yX + ^ 

is also a continuous function of y in L l (B) norm and is a kernel of L B = 
(Id — K B )~ 1 K B - Taking Bj in (1.34) infinitesimally small one concludes that 
for each the distribution of k particles x\ < X2 < ■ ■ ■ < it in 5 has 
a density with respect to the Lebesgue measure. Denoting this density by 
Pk(xi, . . . ,Xk) we obtain 

p k (xi, . . . ,x k ) = det(Id - if B ) • det(L B (a:i,a;j)) 1 < iJ < fc (1-36) 

(It should be noted that (1.36) may be nonnegative even for non-Hermitian 
kernel K, it is easy to see that such K still has nonnegative minors). It 
follows from the definition of /c-point correlation functions that 



Pk(%li • • • , %k) — "T7 / Pk+j( x li • • • j x k, %k+l, • • • , %k+j)dx k+ \ 

j=l J' -> Bj 



1.37) 
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The system of equations can be inversed 



(—iy f 

Pk(x U ...,X k ) = y. 7. / Pk+j(xi, . . . i x ki x k+l, ■ ■ ■ i Xk+j)dXk+l ' ' " dXk+j 

(1.38) 



3=0 



Functions p k (xi, ■ ■ ■ ,x k ) are called Janossy probability densities (see [DVJ], 
p. 122) or exclusion probability densities (see [Ma]). It is easy to check that 



00 I r 

— I Pj(xi, . . . , Xj)dxi . . . dxj = 1 

j=0 3' J B 3 



1.39) 



The r.h.s. of (1.36) still makes sense when \\Kb\\ = ^i(B) = 1 (and therefore 
Pk(xi, . . . , Xk) are properly defined in this case too). Indeed, det(Id — Kb) = 
YYjLi(l — ^j(B)) as a function of Ai has a zero of order 1 at Ai = 1. We 
claim that det(L(xi,Xj))i<ij<k has a pole at Ai = 1 also of order 1. To 

see this we write L = L + L, where L i: j = izr^m ' Vi( x i) ' <Pi( x j), L = 
T.e>2 l-x^B) ■<Pe( x i)-<Pe( x j)- Tlien det(L(x i: Xj))^^ = A k (L(x i: x j ) 1 < ijj < k ), 
and we use the fact that rank (L) = 1. If 1 is a multiple eigenvalue of 
K ■ U B , say \i{B) = \ 2 {B) = ■■■ = \ m (B) = 1 > X m+1 (B) one defines 

— YlT=i i-x^Pb) i Pt( x i) ' Vti. x j) an d proceeds in a similar manner. 

Remark 5. Following Macchi, we call a random point field regular if for 
any Borel B C E satisfying # B < oo (P-a.e.), the generating function Ez #s 
is entire. It follows from our results (see also Theorem 4 below) that any 
determinantal random point field is regular. 

Remark 6. In [Ma] (Theorem 12, p. 113) (see also [DVJ], p. 138) Macchi 
essentially claimed that a necessary and sufficient condition on the integral 
operator K, locally trace class, to define a regular fermion (=determinantal in 
our notations) random point field is < K < 1. As one can see from Theorem 
3 above this condition is sufficient, but not necessary (as we established in 
Theorem 3, the necessary and sufficient condition is < K < 1). For 
completeness it should be noted that Macchi studied the case of continuous 
K(x,y) with TrK < oo. 

Remark 7. Formula (1.36) was established in [Ma], p. 113 (see also [DVJ], 
p. 138 and [TW1], p. 820). 
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We finish §1 with a few more results of general nature about determinantal 
random point fields. 

Theorem 4 

a) The probability of the event that the number of all particles is finite is 
either or 1, depending on whether TiK is finite or infinite. 

b) The number of particles is less or equal to n with probability 1 if and 
only if K is a finite rank operator with rank (K) < n. 

c) The number of particles is n with probability 1 if and only if K is an 
orthogonal projector with rank (K) = n. 

d) For any determinantal random point field with probability 1 no two 
particles coincide. 

e) To obtain results of the theorem for B C E one has to replace K by 
K B . 

Proof of Theorem 4. 

a) One direction is obvious. Indeed, if TrfT = < +oo, then < 
+oo with probability 1. Let us now assume TrK = +oo. Consider a 
monotone absorbing family of compact sets {Bj} ( *L 1 (i.e., Bi C B i+1 
and [JZiBi = E). Then TrK Bj — ► +oo. Fix arbitrary large N. By 

j-*oo 

the construction of {Bj} we have P(#e < N) = lim^oo -P(#b j < N). 
But P(# B3 < N) < 2 N ■ E2-* B i = 2 N ■ det(Id - \ ■ K Bj ) '< 2 N ■ 

b) If rank (K) = n, then writing K(x, y) = Yli=i ^i'fii^'fiiu) (a.e.), and 
p n (xi, ...,x n ) = det(K(xi, Xj))i<ij< n we observe that p m (xi, . . . x m ) = 
(a.e.) for any m > n. Therefore • (#e — 1) • • • • • (#£ — n) — 
f Pn+i( x i, ■ ■ ■ x n+ i)dxi . . . dx n+ i = which implies < n with prob- 
ability 1. 

In the opposite direction, if # E < n (a.e.) we have f Bn +i Pn+i( x i, ■ ■ ■ , 
x n+ i)dxi . . . dx n+ i = for any bounded Borel B C E, therefore Tr(A n+1 
(Kb)) = 0. Since K > we obtain rank (Kb) < n for arbitrary com- 
pact B, which implies rank (K) < n. 
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c) follows from b) and the formula Var(# s ) = Tr(K - K 2 ) = ]\ n i=1 A, • 
(1-Ai). 

d) Let i? n = [—n, n] d . It is enough to show that for any n with probability 
1 no two particles in B n coincide. Let e be arbitrary small. Then 
P{3i ^ j : Xi = Xj G B n } < P{3i ^ j : \x t - Xj\ < e,Xj G B n ,Xj G 

< f Bn (f lx ^ yl<t p 2 (x,y)dx)dy. Since p 2 (x,y) is locally integrable, 
the last integral can be made arbitrary small by letting e — > 0. □ 

The next result gives a criterion for the weak convergence of determinantal 
random point fields. 

Theorem 5. Lei P and P n , n = 1, 2, . . . &e probability measures on (X, B) 
corresponding to the determinantal random point fields defined by the Hermi- 
tian kernels K and K n . Let K n converge to K in the weak operator topology 
and Tt^bKuXb) — ► ^(xbKxb) f or any bounded Borel B C E. Then the 
probability measures P n converge to P weakly on the cylinder sets. 

Proof of Theorem 5. It follows from [SI], Theorem 2.20, p. 40 that the 
assumptions of the theorem imply 

Tr|(X n - K) B \ = \\{K n - ^bW, — . 0. (1.40) 

n— +oo 

As a consequence of (1.40) we have 

Tr • (K n ■ XBl -...-K n . X B m ) — > Tr(K ■ xb 1 • - K • X bJ (1.41) 

n^oo 

for any compact B ± , ... , B m . 

Thus using (1.26), (1.27) one can see that the joint moments of 
with respect to P n converge to the joint moments with respect to P. Since 
the moments of j^s in the case of the determinantal random points define 
the distribution of uniquely one can see (exercise) that P n — P. □ 

The rest of the notes is organized as follows. Section 2 is devoted to the 
various examples of determinantal random point fields arising in Quantum 
Mechanics, Statistical Mechanics, Random Matrix Theory, Representation 
Theory, Probability Theory (Renewal Process, 2D Random Growth Models). 
In §3 we discuss ergodic properties of the translation invariant determinantal 
random point fields. We also point out a special role played by the sine 
kernel K(x,y) = s ™7^~^ . In §4 we discuss the Central Limit Theorem for 
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the counting measure and the Functional Central Limit Theorem for the 
empirical distribution function of spacings. 

It is a great pleasure to thank Ya. Sinai for the encouragement to write 
this paper, B. Simon for the explaination of the result of Lemma 1 , G. 
Olshanski for many valuable remarks, and A. Borodin, B. Khoruzhenko, R. 
Killip, and Yu. Kondratiev for useful conversations. 



2 Examples of Deter minantal Random Point 
Fields 

2.1 Fermion Gas 

Let H = — Jp- + V(x) be a Schrodinger operator with discrete spectrum 
acting on L 2 (E). Let {<fe}eLo be an orthonormal basis of the eigenfunctions, 
Htpt = Xe ■ tpe, Ao < Ai < A 2 < . . . . Consider the n th exterior power of 
H, A n (H) : A n (L 2 (E)) -> A n (L 2 (E)), where A n (L 2 (E)) = A n L 2 (E n )) is 
the space of square-integrable antisymmetric functions of n variables and 
A n (H) = Ya=i( [—£? + V(xi)). In Quantum Mechanics A n (H) describes the 
Fermi gas with n particles. The ground state of the Fermi gas is given by 

if>(xi, ...,x n ) = 

1 n 1 (2 1) 



n - 777 77 vn! 



It could be noted that if)(xi, . . . , x n ) coincides up to a sign e(x 1 , . . . , x n ) with 



the ground state of Yli=i(~ ~£? + V( x i)) acting on S n L 2 (E n ) with the bound- 
ary conditions ip\ Xi=x . = 0. According to the postulate of Quantum Mechan- 
ics the absolute value squared of the ground state defines the probability 
distribution of n particles. We write 

p(x u ...,x n ) = |^(xi,...,x n )| 2 = det ((pi-^Xj)) 

1 nl ~ (2.2) 

• det(^_ 1 (x i )) 1 <. -< n = ^det(X n (^,^))^. .< n , 

where K n (x, y) = Y17=o Vi-iC^OVi-iO/) * s the kernel of the orthogonal projec- 
tor onto the subspace spanned by the first n eigenfunctions of H. We claim 
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that (2.2) defines a determinantal random point field. Indeed, the fc-point 
correlation functions are given by 



(n) ^* f 

Pk V^l) • • • ' X n) = 7 TTT / PnV^l) • • • 5 X n )dXk+\ ■ ■ ■ dx n = 

(n - k)\ J 



(2.3) 



The last equality in (2.3) follows from the general lemma well known in 
Random Matrix Theory. 

Lemma 4 [Me], p. 89). Let (E,d/i) be a measurable space and a kernel 
K : E 2 -> R 1 satisfy 

[ K(x, y) ■ K(y, z)d^y) = K(x, z) (2.4) 

J E 

I K(x, x)d^{x) — const (2.5) 

J E 



Then 



j det(K(x i ,x j )) 1 ^.^ n dii(x n ) = 

J E 

(const - n + 1) • det(K(x i ,x j )) 1< . r 



(2.6) 



We shall consider in more detail two special cases of H . The first case is 
the harmonic oscillator 

a) H = —£s + x 2 , E = R 1 . Then 

Mx) = a 'J'^ exp (t) I? < exp( - l2) » (27) 

are known as Weber-Hermite functions. To pass to the thermodynamic 
limit n^oowe make a proper rescaling 

71 

Xi = ——ryi, i = l,...,n. (2.8) 
(2n) 2 
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Then the Christoffel-Darboux formula and the Plancherel-Rotach asymp- 
totics of the Hermite polynomials ([E]) imply that 

n-l 

K n (x u x 2 ) = y^ J (pt(xi)(pe(x2) = 



(!) 



<Pn(xi) ■ Vn-lM ~ <Pn(x2)<Pn-l(xi) 



X\ - X 2 



has a limit as n — > +00 



K n (x u x 2 ) — > K(y u y 2 ) 



sin7r(yi - 1/2) 



(2.9) 



7r(yi - 1/2) 

The convergence of kernels implies the convergence of /c-point corre- 
lation functions, which in turn implies the weak convergence of the 
distribution 



7T 



(2n 



7T 



7T 



T y n dj/i . . . dy„ 



(2n)2 (2n)s 

to the translation-invariant determinantal random point field with the 

1-2/2) 
-2/2) 



"sine kernel" JT( yi , y 2 ) = 



b) For another example let E = S 1 = {z = e id ,0 < < 2tt},H = -jp. 
Then 

MO) = ' 



Me 



2tt 



10) 



l<j,fc<n 



-det(X n (^,^) 



n! 



l<i J<n 



where 



i^iA) 



1 sin (f • (0 2 - 00) 



(2-11) 



2tt sin (^1) 

After rescaling = y?, i = l,...,n the rescaled correlation functions 
have the same limit as in (2.9), in particular 



2tt 



n— >oo ?7, 



2tt 2tt 



lim — X n — yi, — 1/2 



n 



sin 7T (1/2 - yp 
^2 - yi) 
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For more information we refer the reader to [D1]-[D3], [L4]-[L5], [Sp]. 



2.2 Coulomb Gas at j3 = 2 

Examples a), b) from §2.1 can be reinterpreted as the equilibrium distribution 
of n unit charges confined to the one-dimensional line (ex. 2.1a)) or the unit 
circle (ex. 2.1b)) repelling each other according to the Coulomb law of two- 
dimensional electrostatics. Writing the potential energy as H(zi, . . . , z n ) = 
n log \zi — Zj\ + Y^i=\ V{ z i)-> where V is an external potential, we 



-E 

note that the Boltzmann factor -| exp(— /3H(zi, . . . , z n )), (3 = 2, is exactly 



p n (zi,...z n ) in 2.1a) with V(z) = \z , and p n (9i, . . . , 6 n ) in §2. lb) with 

V{z) (>..;, '"'•./ 1 "• 

The one-component Coulomb gas in two dimensions ( two-dimen- 
sional one-component plasma) was studied in a number of papers including 
[Gin], [Jal], [Ja2], [AL], [DFGIL], [FJ1]. This subject is closely related to 
the theory of non-Hermitian Gaussian random matrices (to be discussed in 
§2. 3d). The two-component two-dimensional Coulomb gas (i.e. a system of 
positively and negatively charged particles) was studied in [Ga], [CJ1]-[CJ3], 
[AF], [FJ2]. Let us start with a neutral system of n positive and n negative 
particles. After denoting the complex coordinates by Uj and Vj,j = 1, . . . , n, 
we write the Boltzmann factor at j3 = 2 as 



exp ( 2 ^ (log \ui-Uj\ + log \vt 



l<i<j<n 



n 



l<i<j<n \ U i U j\ 



Vi - Vi 



Vj\ - 2 log \ui - Vj\) J = 
1 



det 



Ui - v« 



3 / l<i,j<n 



Discretizing the model one allows the positive particles to occupy only the 
sites of the sublattice 7 • Z 2 and the negative particles to occupy only the sites 
of the sublattice 7 • (Z 2 + (|, |)). The grand canonical ensemble is defined 
by the partition function (let 7 = 1) 

1 



\u — v\ 



+ 



2! 



\ + {u 1 )\ + {u 2 )\ + (v 1 )\ + (v 2 ) 



U1,U2,V1,V2 



det 



3' l<i,j<2 



+ 
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where \+(u) = e y( "\A_(-u) = e v< - u \ are fugacities and V is an external 
potential. One can rewrite the last formula as 



where a x , a y , a z are 2x2 Pauli matrices. In particular we see that the grand 
canonical ensemble is a discrete fermion random point field (the appearance 
of matrix- valued kernel reflects the fact that E = Z 2 |J(Z 2 + (|, |))). Passing 
to the continuous limit (7 = 0) one can see that two- and higher order 
correlation functions have a limit, and the limiting kernel K can be expressed 
in terms of the Green function of a differential Dirac operator, namely 

K = I m + ■ h m_ • 



^4 + OyO y + m + V m. 



2 2 

where m + , m_ are rescaled fugacities. In the special case m + = m_ = const 
(i.e., V = 0), K = f^-~'~ ~[~' ^ J can be expressed in terms of modified 



K-+, K 

Bessel function (for the details see e.g., [CJ3]). 



2.3 Random Matrix Models 

a) Unitary Invariant Ensembles of Hermitian Randon Matrices 

The probability distribution in §2. la) (formulas (2.2), (2.7)) allows yet 
another interpretation. It is well known in Random Matrix Theory as the 
distribution of the eigenvalues in the Gaussian Unitary Ensemble (G.U.E.). 
We recall the definition of G.U.E. Consider the space of n x n Hermitian 
matrices {A = (Aj)i<i,j<n, Re(Aj) = BB(A ji ),lm(A ij ) = -Im(A^)}. A 
G.U.E. random matrix is defined by its probability distribution 

P(dA) = const„ • exp(-TrA 2 )dA, (2.12) 

where dA is a flat (Lebesgue) measure, i.e., dA = \\ i< jdRe(A i j) 
dlm^Aij) YYk=i dAkk- The definition of G.U.E. is equivalent to the require- 
ment that {Re(A i j),lm(A i j),l < i < j < n,A kk ,l < k < n} are mutually 
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independent and Re(Aj) ~ N(0, |),Im(Aj) ~ N(0, \),A kk ~ JV(0, §). The 
eigenvalues of a random Hermitian matrix are real random variables. For 
the derivation of their joint distribution we refer the reader to [De], sections 
5.3-5.4 and [Me], chapters 3, 5. It appears that the density of the joint dis- 
tribution with respect to the Lebesgue measure is given exactly by (2.2), 
(2.7). 

We remark that the distribution of a G.U.E. random matrix is invari- 
ant under the unitary transformation A — > UAU~ l ,U G U(n). A natural 
generalization of (2.12) that preserves the unitary invariance is 

P(dA) = const„ • exp(-2 • TrV(A))dA (2.13) 

where V(x) can be, for example, a polynomial of even degree with a positive 
leading coefficients (see [De], section 5). The derivation of the formula for the 
joint distribution of the eigenvalues is very similar to the G.U.E. case. The 
density p n (Xi, . . . , A n ) is given by (2.2), where {<Pe( x )' eV ^}iZo are the first n 
orthonormal polynomials with respect to the weight exp(— 2V(x)). K n (x,y) 
is then again a kernel of a projector and therefore satisfies the conditions of 
Lemma 4. 

b) Random Unitary Matrices 

Let us consider the group of n x n unitary matrices U(n). There exists 
a unique translation invariant probability measure on U(n) (see [We]). It is 
called the Haar measure, we will denote it by /i# aar . The probability density 
of the induced distribution of the eigenvalues is given by 

p n (9 1 ,...,9 n )^(2n)- n -^- J] |e**- e *«| 2 , 

l<k<£<n 

which coincides with (2.10)-(2.11) (see [Me], ch. 9-10, [D1]-[D3]). In the last 
formula we used the notations 

\\ — e , . . . , A n — e . 

If one starts with the probability measure const„ • e~ TrV ^ dfiHaar(U) on the 
unitary group instead of the Haar measure, and replaces the monomials 

-±=e M by ipe(0) ■ e~2 V<y0 \ where {^}"=o are the first n orthonormal poly- 
nomials in e l9 with respect to the weight e~ v( - 9 ^d9, one stills arrives at the 
formula (2.10) for the /c-point correlation functions. 
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c) Random Orthogonal and Symplectic Matrices 

The distribution of the eigenvalues of a random orthogonal or symplectic 
matrix (with respect to the Haar measure) also has a form of a determinantial 
random point field with a fixed number of particles. For the convenience of 
the reader we draw below the chart of the kernels appearing in the ensembles 
of random matrices from the Classical Compact Groups. 





K n {x,y) 


U(n) 




SO(2n) 




SO(2n+ 1) 


1 / sm(n-(x-y)) sin(n- (x+y)) \ . 771 _ rr, 1 
2n [ sin (£_£) sin (£+K) J ' i ' J 


Spin) 


1 f sm(^-(x-y)) si n(^±l-(x+y))\ ^ [n , 
2tt ^ sin(H^) sin(^) ^ ' ^ ~ ^' ^ 


0_(2n + 2) 


the same as for Spin) 



For additional information we refer the reader to [Jol], [DS], [KS], [Sol], 
[So2], [So3]. 

d) Complex Non-Hermitian Gaussian Random Matrices 

In [Gin] Ginibre considered the ensemble of complex non-Hermitian ran- 
dom n x n matrices where all 2n 2 parameters {ReAij, ImAij, 1 < i, j < n} 
are independent Gaussian random variables with zero mean and variance \. 
The joint probability distribution of the matrix elements is then given by the 
formula 

P(dA) = const„ • exp(-Tr(A* • A))dA, 

dA= Y[ dReA jk ■ dlmA jk . ( 2 - 14 ) 

l<j,k<n 

The equivalent definition of (2.14) is that A = A + i ■ A, where A and A are 
two independent G.U.E. matrices. The eigenvalues Ai,...,A n are complex 
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random variables. It was shown that their distribution is given by the deter- 
minantal random point field in M 2 with a fixed number of particles = n) 
and the correlation functions 

Pk\zi, ...,z k ) = de^Knizj^))^.^ (2.15) 

where K n (zi,Z2) = ^exp(— ^ ^-) • Y^i>Zl ^r~- We mention in passing 

that K n (zi,Z2) converges to the kernel 

K( Zl ,zl) = \ exp(-^ - + z x ■ zi) (2.16) 

which defines the limiting random point field. A generalization of (2.14) was 

studied in [Girl], [Gir2], [SCSS], [FKSl], [FKS2]. Let A = A + i • v • A, 

where A and A are, as above, two independent G.U.E. matrices and v is a 
real parameter (it is enough to consider < v < 1). Let us introduce a new 
parameter r = jtti- The distribution of the matrix elements is given by 



P(dA) = const n • exp ( —^—^Tt(A*A - rRe(A 2 )) ) dA (2.17) 



It induces the distribution of the eigenvalues 

n r i n 

p n {zi, , z n ) Y\ dzjdzj = const n • exp - ^ - • ^ 

-i n 



(2-18) 



j<k j=l 



It should be noted that the expression (2.18) also appeared in the papers 
by DiFranceso et al. ([DFGIL]) and Forrester- Jancovici ([FJ1]) as the Boltz- 
mann factor of the two-dimensional one-component plasma. For the calcula- 
tion of the correlation functions we refer the reader to [DFGIL], [FJ1],[FKS1],[FKS2]. 
The crucial role there is played by the orthonormal polynomials in the com- 
plex plane with the weight 



w 2 (z) = exp 



T [z 2 + z 2 ) 



1 -r 2 V 1 1 2 



(2.19) 
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Such orthonormal polynomials can be expressed in terms of the Hermite 
polynomials, 



r 2 



7T2 



(ay. 



(2.20) 



where Y^=o^n{z) ■ = exp(zt — y). We remark that if r = (Ginibre 
case) ipe(z) = j 1 ^ ■ z . Formula for the correlation functions (r arbitrary) 

generalzies (2.15): 



P { k ] =det (K n (zi, Zj)) 1<ij<k , 



n—1 



K n (z 1 ,z 2 ) = w(z 1 )w(z 2 ) ■^ip e (z 1 )ip e (z 2 ) 



(2.21) 



In the limit n — > oo K n (z, z 2 ) converges to 
K(z 1 ,z 2 } = lim X n (zi,^) = 



1 



cxp 



z l\ 



+ 



Z 2 



Z\Z 2 \ 



[2.22) 



tt(1 -r 2 ) 1 - V 2 2 

We reamrk that the last formula differs from (2.16) only by the trivial rescal- 
ing z — > z ■ vl — r 2 . A special regime, called the regime of weak non- 
Hermiticity, was discovered for the model (2.17) by Fyodorov, Khoruzhenko 
and Sommers in [FKS1], [FKS2]. Let 

. . i _i 

Re(zi) = n 2 • x + n ?xi, 

- — - 

Re(z 2 ) = n 2 • x + n 2 x 2) 

Im(zi) = n~^-y 1 , 

lm(z 2 ) = n" 2 • y 2 , 

Assume the parameters x,Xi,x 2 ,yi,y 2 fixed and take the limit n — > oo in 
such a way that lim^oon • (1 — t) = Then 



lim -Kn^zz) = — 

n->oo n 7TQ; 



cxp 



y 2 i + yj , f x (2/1-2/2) 



cr 



S-a I Z I 



^2.23) 
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where 



4 du 



a 2 u 2 



- 2uy 



(2.24) 



(if x > 2 the limit in (2.23) is equal to zero). The formulas (2.23), (2.24) 
define determinantal random point field in R 2 , different from (2.16). 

e) Positive Hermitian Random Matrices 

Following Bronk [Br] we define the Laguerre ensemble of positive Hermi- 
tian n x n matrices. Any positive Hermitian matrix M can be written as 
M = A* A, where A is some complex matrix. The probability distribution of 
a random matrix M is given by 

const„ • exp(-TrAM) • [det (A* A)} a d A, (2.25) 

where dA is defined as in (2.14) and a > — 1 (the values of a of special 
interest are ±|,0). The induced probability distribution of the (positive) 
eigenvalues is given by 

(n \ n 

-J2 Xi ) Tl X 'i II (^-^fdXi-.-Xn (2.26) 
i=l / j=l l<i<j<n 

Employing the associated Laguerre polynominals 

L°L{x) = —e x x- a — (e- x x m+a ), m = 0, 1, . . . 
mV ' n\ dx mK ' 

one can rewrite (2.26) as 

-det( Knfaxj)) (2.27) 

U - \ J l<i,j<n 

where 

n-\ 



K n (x,y) = J2^\x)-^\y), (2-28) 



and {^ a) (x) = (r(ct + 1) • i^))' 2 L e ( x )}e^o is the orthonormal basis with 
respect to the weight e~ x ■ x a on the positive semiaxis. Once again Lemma 4 
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allows us to calculate explicitly fc-point correlation functions and show that 
they are given by determinants of k x k matrices with the kernel (2.28). 

f) Hermitian Matrices Coupled in a Chain 

Let Ai,...,Ap be complex hermitian random n x n matrices with the 
joint probability density 



const n • exp 

+ c x A x A 2 + c 2 A 2 A 3 



-Tr^V^AJ + V 2 (A 2 ) + ■■■ + ^i(Vi) + \v p {A P ) 



c p~iA p -iAp) 



(2.29) 



We denote the eigenvalues of Aj (all real) by Xj = (Xji, ■ ■ ■ , Xj n ),j = 1, 
The induced probability densit of the eigenvalues is then equal to 



,P- 



P n (Ai, . . . , A p ) = const n 



J~J (Ai r — Xi s )(X pr — A ps )] 



Kr<s<n 



■p-1 



Y[det[w k (X kr , X k+ls )] r<s=1 ^ n 



■k=l 



(2.30) 



where 



w k (x,y) = exp ( --V k (x) - ^V k+1 (y) + c k xy 



(2.31) 



Eynard and Mehta [EM] established that the correlation functions of this 
model pfci,...,fc p ( An, • • • , Ai fcl ; . . . ; A pi , . . . , X pkp ) = Y\ P j=i (n-kj)\ I Pn(^u ■ ■ ■ ■> K)' 
Y\ P j=i rir 7 -=fe 7 -+i dXj rj can be written as a k x k determinant with k — ki + 



+ kp 



det [K i:j (X ir , X js )] r=1 



.,ki\s=l,...kj\i,j=l,...,p 



(2.32) 



For the exact formulas for the kernels Kij(x,y) we refer the reader to [EM] 
(see also [AM]). We remark that (2.32) defines a determinantal random point 
field with one-particle space E being the union of p copies of R 1 . 

g) Universality in Random Matrix Models. Airy, Bessel and sine Random 
Point Fields. 
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We start with a general class of kernels of the form 



(2 .33) 

x-y 

where 

m(x)ip'(x) = A(x)ip(x) + B(x)i)(x) 
m(x)ip'(x) — — C(x)ip(x) — A(x)ip(x) 

and m(x),A(x),B(x),C(x) are polynomials. It was shown by Tracy and 
Widom ([TW2]) that Fredholm determinants of integral operators with ker- 
nels (2.33)-(2.34) restricted to a finite union of intervals satisfy certain par- 
tial differential equations. Airy, Bessel and sine kernels are the special cases 
of (2.33), (2.34). To define sine kernel we set (p(x) = ^sin(7rx), ip(x) = 
(p'(x) (m(x) = 1, A(x) = 0, B(x) = 1, C(x) = ir 2 ). For the Airy kernel 
ip(x) = Ai(x), ip(x) = (p f (x) (m(x) = 1, A(x) = 0, B(x) = 1, C(x) = —x). 
For the Bessel kernel (fi(x) = J a (\/x), ip(x) = xip'(x) [m(x) = x, A(x) = 
0, B(x) = 1, C(x) = \{x — a 2 )). Here Ai(x) is the Airy function and J a (x) is 
the Bessel function of order a (see [E]). Writing down these kernels explicitly 
we have (see [TW2], [TW3], [TW4]) 

K,„4 X ,y)= 3m ; {x -«K (2.35) 



K Airy {x,y) = 



Jo 



x-y 



(2.36) 
Ai(x + t)Ai(y + t)dt 



K B e SSe l{X,y) = — r 

{X ~ V) ^ (2.37) 

y/x- Jg+ijVx) ■ J a {y/y) - J a {y/x)y/y- Jg+iiy/y) 

2- (x-y) 

As we already mentioned, sine kernel appears as a scaling limit in the bulk of 
the spectrum in G.U.E. ([Me], Chapter 5). In its turn, Airy kernel appears as 
a scaling limit at the edge of the spectrum in G.U.E. and at the (soft) right 
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edge of the spectrum in the Laguerre ensemble, while Bessel kernel appears 
as a scaling limit at the (hard) left edge in the Laguerre ensemble ( [F] , [TW3] , 
[TW4]). Universality conjecture in Random Matrix Theory asserts that such 
limits should be universal for a wide class of Hermitian random matrices. 
Recently this conjecture was proven for unitary invariant ensembles (2.13) in 
the bulk of the spectrum ([PS], [BI], [DKMLVZ], [De]) and for some classes 
of Wigner matrices in the bulk of the spectrum [ Jo4] and at the edge [So4] . 

In the next subsection we completely characterize determinantal random 
point fields in M 1 (Z 1 ) with independent identically distributed spacings. 

2.4 Determinantal random point fields with i.i.d. spac- 
ings. Renewal processes 

We start with some basic facts from the theory of renewal processes (see e.g., 
[Fe], [DVJ]). Let {t^}^ be independent identically distributed non- negative 
random variables and To another non-negative random variable independent 
from {r fe }^ 1 (in general the distribution of r will be different). We define 

k 

x k = Ys T r ( 2 - 38 ) 

3=0 

This gives us a random configuration {xk}kL in R+. In probability theory a 
random sequence {xk}^L Q is known as delayed renewal process. We assume 
that the distribution of random variables Tk,k > 1 has density f(x), called 
interval distribution density, and a finite mathematical expectation Et\ = 
Jo°° x f( x )dx. The renewal density is defined than as 

<x) = f k *(x) = f(x) + f f{x - y)f{y)dy + 

fc=1 (2.39) 

/ / f{x-y\- y2)f(yi)f(y2)dyidy 2 + ... 

Jo Jo 

One can express higher order correlation functions of the renewal process 
through its one-point correlation function and the renewal density. Indeed 
(see [DVJ], p. 136) for t± < ti < • ■ ■ < tk and k > 1 the following formula 
takes place 

p fc (fi, . . . , f fc ) = pi(fi) • u(t 2 - ti) • u(t 3 - t 2 ) ■ . . ■ ■ u(t k - t fc _i) (2.40) 
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It follows immediately from the above definitions that a random point field 
in M.+ has i.i.d. nearest spacings iff it is a renewal process (2.38). To make 
this process translation-invariant the probability density of To must be given 
by 

i r+oo 

/ f(t)dt ([DVJ], p. 72, [Fe], section XI.3) (2.41) 

En J x 

Then one-point correlation function is identically constant, pi(x) = p > 0, so 
(2.40) implies that the distribution of the process is uniquely defined by the 
renewal density (in particular one can obtain p from u(x) since p = (Eti)^ 1 
and the Laplace transforms of / and u are simply related). Macchi ([Ma]) 
considered a special class of translation-invariant renewal processes with the 
interval distribution density 

f(x) = 2p(l-2pa)~h~^ -smli^l - 2pa)^ ■ {^fj, (2-42) 

where 

2pa < 1, p > 0, a > 0, (2.43) 

and showed that it is a determinantal random point field with the kernel 

K(x,y) = p-exp{-\x-y\/a) (2.44) 

(restrictions (2.43) are exactly < K < Id). 

In the next theorem we classify all delayed renewal processes that are also 
determinantal random point fields in . 

Theorem 6. Determinantal random point field in with Hermitian ker- 
nel has i.i.d. spacings if and only if its kernel satisfies the following two 
conditions in addition to < K < Id locally trace class : 

a) for almost all x\ < x 2 < x 3 

K(x 1 ,x 2 ) ■ K(x 2 , x 3 ) = K(x 1 ,x 3 ) ■ K(x 2 , x 2 ), (2.45) 

b) for almost all x\ < x 2 the function 

K(X2 ' X2) K^) (2 ' 46) 
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depends only on the difference x 2 — X\. If a determinantal random point field 
is both translation-invariant and with i.i.d. spacings, it is given by (2.42)- 

(2.44)- 

Remark 8. Of course a translation-invariant d.r.p.f. in M.+ can be extended 
in a unique way to the translation-invariant d.r.p.f. in R 1 . 

Proof of Theorem 6. First we prove the "only if part of the theorem. 
Suppose that a determinantal random point field with a kernel K(x, y) is 
also a delayed renewal process. From (2.40), k = 2, 3, we obtain the formula 
for the renewal density 

u(y -x)= K(y, y) - V > x, (2.47) 

and the expression for p 3 (xi,x 2 ,x 3 ),x 1 < x 2 < x 3 : 

p 3 (x 1 ,x 2 , x 3 ) = K(x 1 ,x 1 ) ■ u(x 2 - xi) ■ u(x 3 - x 2 ) 

K(x 1 ,x 2 ) ■ K(x 2 ,xi) 



= K(x 1 ,x 1 ) ■ [ K(x 2 ,x 2 ) 
K(x 3 ,x 3 ) 



K( Xl , Xl ) J (2.48) 
K(x 2 ,x 3 ) ■ K(x 3 ,x 2 ) 



K(x 2 ,x 2 ) 

Since with probability 1 there are no particles outside A = {x : K(x, x) > 0}, 
we can always consider random point field restricted to A. 
Comparing 

p 3 (x 1 ,x 2 , x 3 ) = det(K(xi, x_,-))i<ij<3 (2.49) 
with (2.48), we have 

K(x!,x 2 ) ■ K(x 2 , xi) ■ K(x 2 , x 3 ) ■ K(x 3 , x 2 ) ■ ——^ = 

K{x 2 ,x 2 ) 

- K(xi,x 3 ) ■ K(x 3 , X!) ■ K(x 2 , x 2 ) + K(x!,x 2 ) ■ K(x 2 , x 3 ) 
■ K(x 3 , xi) + K(xi,x 3 ) ■ K(x 3 , x 2 ) ■ K(x 2 , xi) 

which is equivalent to 

1 

K(x 2 ,x 2 ) 

%1 „|, ( , I , lI1 )- %I ,), ( , IllI ,)H 



K(x!,x 2 ) ■ K(x 2 , x 3 ) - K(x 2 , x 2 ) ■ K(x!,x 3 ] 
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The third factor in the last equality is a complete conjugate of the second 
factor, and we obtain (2.45). Condition b) of the theorem has been already 
established in (2.47). For the translation-invariant d.r.p.f. the kernel K(x,y) 
depends only on the difference, therefore K(x, y) — p- e~^ x ~ y ^ a ■ e l ^ x ~ y \ and 
the unitary equivalent kernel e~ l(3x K(x,y) ■ e l/3y coincides with (2.44). Now 
we turn to the proof of the "if part of the theorem. Once we are given the 
kernel satisfying (2.45) and (2.46) the candidate for the renewal density must 
obey u(x 2 — x±) = K(x 2 ,x 2 ) — K ( X1 ^' K ( X2 ' X1 ^ f or almost all x\ < x 2 . Let 
x\ < x 2 < ■ ■ ■ < Xk- Our goal is to deduce the algebraic identity 

fe-i 



det ( K(xi, Xi) I = K(x 1 ,x 1 ) ■ JJ ( K(x i+1 ,x i+1 ) 

^ ' l<i,j<k 

K(xi,x i+1 ) ■ K(x 

i+lj x i 



(2.50) 



from the basic identities between the commuting variables K(xi, Xj), K(x i: Xj) 
satisfying K(xi,Xj) ■ K(xj,xt) = K(xi,xt) ■ K(xj,Xj),l < i < j < t < 
k, K(xi,Xj) = K(xj,Xi). Let us introduce a(x) = K(x,x) ■ K(0, x)' 1 , b(x) = 
K(0,x)-\ Then for i <j 

K(xi,Xj) = a(xi) ■ b(x j )~ 1 , 



K(xj,Xi) = a(xi) ■ b(xj) 
This allows us to write the determinant as 



a(xi) ■ b(xi) 1 , 
— i 



a(xi) • b(x 2 ) 



a(xi) ■ b(x 2 ) 
a(x 2 ) ■ b(x 2 ) 



a(xi) ■ b(x n ) 1 
,a(x 2 ) ■ b(x n )~ 1 



a(xi) • b(x n ) , a(x 2 ) ■ b(x n ) 



, a(x n ) ■ b{x n ) 



-i 



n-l 



a{x 1 )-b(x 1 ) 1 • Y[( b(x i+1 ) 1 -b(xi +1 ) 



— i 



(2.51) 



(a(x i+1 ) ■ b(x i+1 ) - a(xi) ■ bfa)) 



which is exactly the r.h.s. of (2.50). Once we established Pk(xi, ■ ■ ■ ,Xk) = 
Pi( x i) ' ni^i 1 u ( x i+i ~ x i)i x i — x 2 < • ' ' < x ki the rest of the proof is quite 
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easy. Let p k (xi, ■ ■ ■ ,Xk) be Janossy densities, i.e. the probability density of 
the event to have particles at x±, . . . , x k , and no particles in between. We re- 
call that Pk(xi, . . . , x k ) = Y^jLi —}¥-J Pk+j( x i, ■■■,x k ; y k+1 , . . . , y k+j )dy k+l . . . dy k+j , 
where the integration in the jth term is over (xi,x k ) x • • • x (xi,x k ). We 

< — j times > 

claim that 

k-l 

p k (x u ...,x k ) = pi(xx) • JJ/(x i+ i -Xi) (2.52) 

i=i 

where the interval distribution density / and the renewal distribution density 
u are related via the convolution equation 

u = f + u*f (2.53) 

Theorem 6 is proven. □ 

Remark 9. The analogue of Theorem 6 is valid in the discrete case and 
the proof is the same. One has to replace (2.42) by the solution of the 
discrete convolution equation (2.53) with u{n) = 1 — p ■ e~ 2/3n , K(ni, n 2 ) = 
p ■ e-fl" 1 "" 2 !, < p < 1, > 0, so that 

^ H - n) - (p~ 213 - n) ■ p u 

m = Y. /(»)»'"' - (2 _ p) i (2 ^- P+ i) Z. + e -v < 2 - 54 ' 

Remark 10. One can consider a generalization of Theorem 6 to the case 
when the multiplicative identity (2.45) still holds, but the renewal density 

v(r T )- K(r r ) _ K(x u x 2 ) ■ K(x 2 , Xl ) 

U{Xi,X 2 ) — J\{X2,X2) 777 \ 

K(x 1 ,xi) 

no longer depends only on the difference of X\ and x 2 - Such processes have 
independent but not necessarily identically distributed spacings since the 
spacings distribution f(x, y)dy depends on the position x of the left particle. 
Thus 

f(x!,yi) ■ f(y 1 ,x 2 )dy 1 
rx2 ryi (2.55) 
+ / / f(xi,yi) • f(yi,V2) ■ f(V2,x 2 )dy 1 dy 2 + ... 

J XI J X1 
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where f(x,y) is a one-parameter family of probability densities, such that 



supp f(x,-) C [x,+oo], f > 0, J f(x,y)dy = 1. 
We recall the inversion formula for the equation (2.55): 
f(xi, x 2 ) = u(xi, x 2 ) - / u(x 1 ,y 1 )u(y 1 ,x 2 )dy 1 



<X\ 



PX2 PV2 

+ / u(xi,yi) ■u(y 1 ,y 2 ) ■u(y 2 ,x 2 )dy 1 dy 2 - ... 

J X-\ J XI 



(2.56) 



Writing K(x,y) = a(x)b(y) \x < y, where a(x) = ^0j,b(y) = 

and u(x,y) = j^yp • { a {y)b{y) — a(x)b(x)), one can in principle characterize 
through (2.56) the class of corresponding interval densities u(x,y). 



2.5 Plancherel Measure on Partitions and its Generali- 
zations— z-Measures and Schur Measures 

By a partition of n — 1, 2, . . . we understand a collection of non- negative 
integers A = (Ai, . . . , A m ) such that Ai + - • -+A m = n and Ai > A2 > • • • > A m . 
We shall use a notation Par(n) for the set of all partitions of n. For the 
basic facts about partitions we refer the reader to [St], [Fu], [Mac], [Sa]. In 
particular recall that each partition A of n (denoted A h n) can be identified 
with a Young diagram with |A| — n boxes. A partition A' corresponds to 
the transposed diagram. Let d be the number of diagonal boxes in A (i.e., 
the number of diagonal boxes in the Young diagram corresponding to A). 
We define the Frobenius coordinates of A as (pi, . . . ,Pd\qi, ■ ■ ■ , Qk), where 
Pj = Xj — j, qj = X'j — j, j = 1, • • • ,d. The importance of partitions in 
Representation Theory can be most easily understood from the fact that 
the elements of Par(n) label the irreducible representations of the symmetric 
group S n (see e.g., [Sa], [Fu]). The Plancherel measure M n on the set Par(n) 
of all partitions of n is given by 

M„(A) = (2.57) 

where dim A is the dimension of the corresponding representation of S n . The 
dimension dim A can be expressed in terms of the Frobenius coordinates via 
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a determinantal formula 



dim A 



= det 



1 



(2.58) 



(Pi + qj + 1) • Pi ] - ■ ft! 



l<i,j<d 



where |A| = n ([01], Proposition 2.6, formula (2.7)). Let Par = |_I^L Par(ra). 
Consider the measure M 9 on Par, which in analogy with statistical mechanics 
can be called the grand canonical ensemble: 



M 9 is also called the poissonization of the measures M n . It follows from (2.59) 
that | A | is distributed by the Poisson law with the mean 9, and M 9 { | |A| = 
n) = M n . In the Frobenius coordinates measures M 9 ,M n can be viewed 
as random point fields on the lattice Z 1 . Recently Borodin, Okounkov and 
Olshanski ([BOO]) and, independently, Johansson [Jo3]) proved that M 9 is a 
determinantal random point field (to be exact in [Jo3] only the restriction of 
M e to the first half of the Frobenius coordinates (pi, . . . ,Pd(\)) was studied, 
and as a result, only the part of (2.60) corresponding to xy > was obtained). 
To formulate the results of [BOO], [Jo3] we define the modified Frobenius 
coordinates of A by 



M e (X)=e- e --M n (\), 



(2.59) 



if A E Par(n), n = 0, 1, 2, ... , < 9 < oo. 



Fr(A) := { Pl + -,..., p d + -,- qi - -q d - -}. 



Let p 9 k (xi, . . . ,Xk) be the fc-point correlation function of M 9 in the modified 
Frobenius coordinates, where 





(2.60) 
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where J x {-) is the Bessel function of order x. We note that the kernel K(x, y) 
is not Hermitian symmetric, however the restriction of this kernel to the 
positive and negative semi-axis is Hermitian. (2.60) can be seen as a limiting 
case of a more general theorem obtained by Borodin and Olshanski for the 
so-called z-measures (see Theorem 3.3 of [BOl], also [B02], [B03], [KOV] 
and references therein) . Let z, z' be complex numbers such that either 

z =zeC\Z 

or (2.61) 
[z] < mm(z, z') < max(z, z') < [z] + 1, 

where z, z' real and [ ] denotes the integer part. Let (x)j — x ■ (x + 1) • . . . • 
(x+j — 1), (x)q = 1. Below we introduce a 2-parametric family of probability 
measures M^, on Par(n). These measures take their origin in harmonic 
analysis on the infinite symmetric group ([KOV], [01]). By definition 



(z ■ z') d[ X) 
(z ■ z% 



= ^j^r Eft* + !)« • V + !)« 

*=i (2.62) 

i t \ i i ,\ dim 2 A 
<-z + l) qi .(-z' + l) qi -——- 



The conditions on z, z' stated above are equivalent to the requirement that 
(z)j ■ (z')j and (—z)j ■ (-z')j are positive for any j = 1, 2, . . . We note that 
M^l converges to the Plancherel measure M n if z, z' — > oo. The measure 

M^, is called the n-th level z-measure. Consider now the negative binomial 
distribution on the non-negative integers 

(l-er'-^^T' n = 0,l,... 

where £ is an additional parameter, < £ < 1. The corresponding mixture 
of the n-level z-measures defines measure M ZyZ >^ on Par. We remark that 
M z ,z'£ degenerates into M e if z, z' — > oo, £ — > in such a way that zz'^ — > ^. 
It was shown in [BOl] that in the modified Frobenius coordinates M z>z >^ is 
a determinantal random point field on Z 1 + \ . The corresponding kernel can 
be expressed in terms of the Gauss hypergeometric function and is called 
the hypergeometric kernel. It appears that a number of familiar kernels can 
be obtained in terms of the hypergeometric kernel, in particular Hermite 
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kernel ((2.2), (2.3), (2.7)), Laguerre kernel ((2.2), (2.28)), Meixner kernel 
((2.67) below), Charlier kernel. For the hierarchy of the degenerations of the 
hypergeometric kernel we refer the reader to [B02] §9. Recently Okounkov 
[Okl] showed that the measures M ZjZ >£ are the special case of an infinite 
parameter family of probability measures on Par, called the Schur measures, 
and defined as 

M(\) = -s x (x)-s x (y), (2.63) 

z 

where s x are the Schur functions (for the definition of the Schur functions 
see [St] or [Mac]), x — (x±, X2, ■ ■ • ) and y — (y±, y 2 , . . . ) are parameters such 
that 

z = s *( x ) ■ s ^(y) = Hi 1 - x ^y l ( 2 - 64 ) 

AS Par i,j 

is finite and {x i } c *^ l = Measures M Z;Z /^ formally correspond to 

x™ = £~ ■ z, Yli^i UT = £ T ' z ' i m = 1) 2 . . . .To be precise one should 
consider the Newton power sums as real parameters and express the Schur 
functions as polynomials in the power sums. By now the reader probably 
would not be very surprised to learn that the Schur measures also can be 
considered as determinantal random point fields ([Okl], Theorems 1,2)! 

2.6 Two-Dimensional Random Growth Model 

As our last example we consider the following two-dimensional random growth 
model ([Jo2]). Let {a>ij}ij>i be a family of independent identically dis- 
tributed random variables with a geometric law 

p( aij = k)=p-q k , k = 0,1,2,... (2.65) 

where < q < l,p — 1 — q. One may think about (2.65) as the distribution 
of the first success time in a series of Bernoulli trials. We define 

G(M, N) = max V a^, (2.66) 

(i,j)£ir 

where the maximum in (2.66) is considered over all up/right paths it from 
(1,1) to (M, N), in other words over vr = {(ii, ji) = (1, 1), (i 2 , J2), J3), 
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(im+n-iJm+n-i) = (M,N)}, such that (i k +i, jk+i)-(h, jk) G {(0, 1), (1, 0)}. 
We mention in passing that distribution of random variables {G(M, N)} can 
be interpreted in terms of randomly growing Young diagrams and totally 
asymmetric exclusion process with discrete time (for the details see [Jo2]). 
Without loss of generality we may assume M > N > 1. To state explicitly 
the connection to the determinantal random point fields we introduce the dis- 
crete weight w 9 K (x) = ( X+I ^~ 1 ) -q x ,K = M — N + 1, on non- negative integers 
x = 0, 1, 2 . . . The normalized orthogonal polynomials {M„(a;)} n >o with re- 
spect to the weight w q K are proportional to the classical Meixner polynomials 
([Ch]). The kernel 



3=0 

satisfies the conditions of Lemma 4 with respect to the counting measure on 
non-negative integers. Therefore 



defines a discrete determinantal random point field. It was shown by Johans- 
son that the distribution of the random variable G(M, N) coincides with the 
distribution of the right-most particle in (2.68). After appropriate rescaling 
in the limit TV— >oo,M— >oo,^— > const, this distribution converges to 
the distribution of the right-most particle in the Airy random point field 
(2.36). Additional information on the subject of the last two subsections can 



be found in the recent papers/preprints [AD], [BDJ1], [BDJ2], [BR1], [BR2], 
[BR3], [Bor], [ITW], [Ku], [Ok2], [PS1], [PS2], [TW5], [TW6] . 



3 Translation Invariant Determinantal Ran- 
dom Point Fields 



As before (X, B, P) denotes a random point field with a one-particle space 
E, hence X is a space of locally finite configurations of particles in E, B is a 
Borel a-algebra of measurarble subsets of X and P is a probability measure 
on (X, B). Throughout this section we always assume E = R d or Z d . We 
define a continuous action {T l } t &E of £ on X in a natural way: 



JV-l 




(2.67) 





T* : X - X, (T'O* = (Oi + t. 



(3.1) 
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Definition 5. Random point field (X, B, P) is called translation invariant if 
for any A G B, any t G E 



P(T-*A) = P{A). 



The translation invariance of a random point field implies the invariance 
of fc-point correlation functions: 

p k {x x + t, . . . , x k + t) = p k (x 1: . . . , x k ), a.e. 

k = 1,2,...,* G ^ ' ' 

Conversely, if {p k } are invariant under {T*}, then there exists a correspond- 
ing random point field which is translation invariant ([L3]). In particular, 
if the translation invariant correlation functions define P uniquely then the 
random point field is translation invariant. In the case of a determinan- 
tal random point field this implies the following criterion: a determinantal 
random point field is translation invariant if and only if the kernel K is 
translation invariant, i.e., K(x,y) = K(x — y,0) =: K(x — y). In this section 
we restrict our attention to the translation invariant determinantal random 
point fields. We are interested in the ergodic properties of the dynamical 
system (X, B, P, {T*}). For the convenience of the reader recall some basic 
definitions of Ergodic Theory ([CFS]). 

- A dynamical system is said to be ergodic if the measure P(A) of any 
invariant set A equals or 1. 

- A dynamical system has the mixing property of multiplicity r > 1 if 
for any functions / , /i, . . . , f r G L r+1 (X, B, P) we have 

lim / /o(0/i(^0 : • • • : f r (T h+ - +t ^)dF = f[ I ftftdP 

(3.3) 

- A dynamical system has an absolute continuous spectra if for any / G 
L 2 (X, B, P) orthogonal to constants 

J f(OW¥)dP = J e^h f {\)d\- (3.4) 
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where the integration at the r.h.s. of (3.4) is over M. d in the continu- 
ous case and over [0,27r] d in the discrete case, and hf(X)d\ is a finite 
measure absolutely continuous with respect to the Lebesgue measure. 
One can interpret (3.4) in the following way. We define a <i-parameter 
group of unitary operators {C/*} t6 E on L 2 (X,B,P) as 

(u t f)(0 = f(T t O- 

Usually such family of unitary operators is called adjoint to the dynam- 
ical system. It is easy to see that {U 1 } commute. Since L 2 (X, B, P) 
is separable and {U 1 ^, ip) is a measurable function of t for any tp,ip e 
L 2 (X,B,P) one can apply the von Neumann theorem ([RS], vol. 1, 
Theorem VIII. 9) to conclude that U l is strongly continuous. In the case 
E = M, d one has hf(X)dX = d(f, Q\f), where dQ\ is a projection- valued 

measure, Q\ = <5(-oo,Ai)x-x(-oo,A d ) = ri f=i X(-oo,A,)(^4j), {A?}j=i are 
the generators of the one-parameter groups £/(°>--->*:i> '--- ) an d X(-oo,t) is 
the indicator of (— oo, t) ([RS], vol. I, Theorem VIII. 12). In the discrete 
case E = Z d dQ\ is a projection- valued measure on a rf-dimensional 
torus, 

d 

Q[l > e«i]x...x[l^] = tlx^V&Uj = 

Theorem 7. Let (X, B, P) be a translation invariant determinantal random 
point field. Then the dynamical system (X, B, P, {T*}) is ergodic, has the 
mixing property of any multiplicity and its spectra is absolutely continuous. 

Remark 11. Recall that the absolute continuity of the spectra implies the 
mixing property of multiplicity 1, which in turn implies ergodicity ([CFS]). 

Proof of Theorem 7. We note that the linear combinations of 

N 

i=i (3.5) 

n>i, s g (0 = J2g(xi),9i e c °°(r*), j = i,...n 

i 

are dense in L 2 (X,B,P). Therefore it is enough to establish (3.3), (3.4) 
for the functions of such form. We start with the lemma calculating the 
mathematical expectation of (3.5). 
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Lemma 5. 

a) 



N 



N 



Mto<0 = E E II E E E 



m=l over partitions 

\J?=iC e ={l,..,N} 



= 1 = 1 o?;er partitions 



(3.6) 



w/iere ^(x) 



/ n^w^O ' _ ^)cfei • • • dx ht 



b)EU 



J'=l 



•+JVr+l O /V\ TT»*+1 



1 - rCi (EnS5::S_ 1+ i S ai (tj) = (similar 



ex- 



pression to (3.6), with the only difference that partitions 

m f r+1 



(3.7) 



e=i 



satisfy (*)), where 

(*) There exists at least one element Cg of the partition such that the 
intersections ofCg with at least two of the following sets {1, . . . N±}, . . . , {iV"i + 
• • - + 7V s _i + l, . . . , N ± + - ■ - + N S }, . . . , {N ± + - ■ - + N r + 1, . . . , JVi + - • - + N r+1 } 
are non-empty. 

Proof of Lemma 5. The proof of part a) is rather straightforward and 
quite similar to the one given at the beginning of §2 in [So3] (see formulas 
(2.1)- (2. 7) from the reference). The proof of part b) follows from a). □ 

To derive the mixing property (3.3) we replace gj(-) for Ni H h N 8 _i + 

1 < j < N 1 + --- + N s ,s= l,...,r + l, in (3.7) by &■(• + *i + • • • + t 8 _i). Fix 
a partition \J^ =1 @t = {1, 2, . . . , A^x + • • • + N r+1 }. Since {gj} are bounded 
functions with compact support, each of m factors at the r.h.s. of (3.7) 
is bounded. We claim that the I th factor (corresponding to Cg, where £ 
is the same index as in (*)) goes to zero. To see this we fix an arbitrary 
partition of Cg, -B« = Cg. By assumption, Cg contains indices 1 < 
u < v < N 1 + ■ ■ ■ + N r+ i, such that u and v belong to different subsets 
{1,...,N 1 },...{N 1 + --- + N s _ 1 + 1,...,N 1 + --- + N S },...,{N 1 + --- + N r + 
1, . . . , iVi H h N r+1 }. We claim that 



/ II^MO^) ' K ( Xi + 1 - x i) dx l ■ ■ ■ dx ke 
J i=l 



(3.8) 
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goes to zero as min{t s , 1 < s < r} — > oo. Indeed if min{t s , 1 < s < r} is 
sufficiently large, the indices u, v belong to different B^s or the corresponding 
gB u is zero (the supports of the factors in gs H will not intersect). Once u 
and v belong to different B^s the argument in K(xi + i — Xi) for some i is 
greater than mm{t S) 1 < s < r}. Since the Fourier transform of K(x), K(t) = 
J e %xi K{x)dx is a non-negative integrable function (bounded from above by 
f), applying the Riemann-Lebesgue lemma we obtain that K(x i+ i —Xj) goes 
to zero. The other terms in (3.8) are bounded and the integration is over a 
bounded set, therefore (3.8) goes to zero and the proof of the mixing property 
follows. 

To establish the absolute continuity of the spectrum we apply (3.7) when 
r = 2, N 1 = N 2 = N, g N+j (x) = 9j (x + t),j = 1,...,N, /(£) = Uf=i S gj (£), 

fiT'O = nL %(^0 = U-= N+1 %(0- We have 



2N 



e(/(o - m ■ - e/) = J2 E HE 



m=l over partitions £=1 

Ufcl^={l,..,2JV} 



over partitions <J^S & 

U k iU Bti=Ct 



{-iy 



r ke . 

/ II 9B la{i) (x^ ■ K(x i+1 - Xi)dx-L . . . dx ke } 
J i=i 



(3.9) 



where we assume that = x\ in the integral, and the sume in J2* i s over 
partitions {Ci, . . . , C m } such that for at least one element Ce of the partition 
both Qn{l,2, ... ,A^} and Qfl{JV + l, . . . , 2iV} are non-empty (we denoted 
above this property by (*)). The terms in the product YYeLi corresponding 
to those £ that do not satisfy (*) are constants as functions of t. Fix now £ 
satisfying (*). We claim that 



IT* 



MO 



1 \ h k e 

I) n 



[xi) ■ K(x i+ i -Xi)dxi... dx ke 



9B ea(i) (Vi+i ~ Vi) ■ K(y i+ i)dy 1 . . . dy k( 



(3.10) 



can be written as J e l ^' x ^h(\)d\, where h(X) is an integrable function. The 
check is rather straightforward and we leave the details to the reader. We 
infer that (3.9) is a linear combination of the products of the Fourier trans- 
forms of integrable functions. Since the product of the Fourier transforms is 
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the Fourier transform of the convolution the proof of the absolute continuity 
of the spectrum follows. Theorem 7 is proven. □ 

One can without difficulty calculate the spectral density of the centralized 
linear statistics 

S g (Z)--ES g = Y^gfa) -Ej2g(xi). 

i i 

Namely 

E(S g - ES^S^ -ES^=J ■ (K(0) - \K\ 2 (X)) 



(3.11) 



-\g(\)\ 2 d\, and h Sa (X) = (K(Q) - \K\ 2 (X)) • -\g(X)\'< 



We conclude that 

H(d\) = (K(0) - \K\ 2 (X))dX (3.12) 

is the spectral measure of the restriction of {U 1 } to the subspace of the 
centralized linear statistics. Since < K(X) < 1, K(0) — J K(X)dX, we 
see that 

< ^ = K(0) - \K\ 2 (X) = K(Q) -±J K(y)K(y - X)dy < K(0). 

We note that g > for A ^ 0, and g(0) = if and only if K(X) is an 
indicator. In particular the spectral measure \i is equivalent to the Lebesgue 
measure. 

Before we formulate the next lemma recall that by #[_L,L] d (0 we denote 
the number of particles in [— L,L] d . 

Lemma 6. 

Var(# hLiL]d )=Vol([-L,L] d )- ^(0) + o(l)) as L - oo. (3.13) 



Proof of Lemma 6. The probabilitists are well familiar with the ana- 
logue of this result in the Theory of Random Processes: let {f] n } be L 2 - 
stationary random sequence and h(X) its spectral density, Er] n rf^ = b(n — 
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m ) = h iT e iX < n - m ^h{X)dX, thenV&Y(ri n + --- + ri n ) = (h(0) + 5(1)) ■ n ([IL], 
section XVIII. 2). To prove the lemma we write 



Var(# [ _ i>L]d ) = / / p2(x,y) - pi(x)p!(y)dxdy 

J{~L,L]*J{-L,L]* 

+ / pi(x)dx — — / / \K\ 2 (x — y)dxdy 

J[-L,L] d J[-L,L] d J[-LM d 

+ AT(0)Vol([-L, L] d ) = (K(0) - [ \K\ 2 (x)dx + o(l)) • Vol([-L, L] d ) = 
(K(0)-\K\i(0) + o(l)).Vo\([-L,L] d ). 

□ 

The sub leading terms in (3.13) also depend on the behavior of ^ near 

the origin. For example, let K(X) be an indicator, K(\) = xb(X),B C W 1 . 
As we have seen above this is equivalent to (0) = 0. For simplicity we will 
assume d—1. If B is a union of m disjoint intervals 

^(A) = ^(0)-^ J K(y) ■ K(y - X)dy = 

— • [length(B) - length^ n {B + A))] = ( 3 - 14 ) 
2% 

^.|A|- (l + o(l)), A-0, 

and after more careful evaluation of the asymptotics of J^ L J^ L \K\ 2 
y )dxdy =±-JZo ■ C sin i L ' X) ) 2 d\ we arrive at 



x 



m 



Var(#[_ L;L]d ) = — logL • (1 + o(l)) (3.15) 



7T 



Choosing m — 1, -R'(A) = X[_ 7r 7r ](A) one obtains the sine kernel K(x — y) — 
• A special role played by the sine kernel can be highlighted by the 
fact that 4j logL rate of the growth for Var(4f[-L,L]) is the slowest among all 
translation-invariant kernels K(x — y) corresponding to projectors, K = xb, 
for which inf(B), sup(B) are the density points of B (if K is not a projector 
it follows from Lemma 6 that the rate of the growth of the variance is linear) . 
As an example consider B = |_L>i[ n ' n + ^1>7 > 1> than one has ^ ~ 



a- 1 



t and Var(#[_ LL] ) ~ Lt. More generally, ^ ~ |A| a ,0 < ct < 1, implies 
Var(# hLiL] ) ~ L 1 -*.' 
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4 Central Limit Theorem for Counting Func- 
tion and Empirical Distribution Function 
of Spacings 

In [CL] Costin and Lebowitz proved the Central Limit Theorem for 
in the case of the sine kernel. The article also contains a remark on p. 71, 
due to Widom, that the result holds for a larger class of Random Matrix 
models. In its general form this theorem appeared in [So2] . 

Theorem 8. Let E be as in (1.1), {0 < K t < 1} a family of locally trace class 
operators in L 2 (E), {(X, B, P t )} a family of the corresponding determinantal 
random point fields in E, and {I t } a family of measurable subsets in E such 
that 

Var t # /t = Tr(K t ■ X i t - (K t ■ X i t f) - oo as t - oo. (4.1) 

Then the distribution of the normalized number of particles in I t (with respect 
to P t ) converges to the normal law, i.e., 

Remark 12. It was shown in [So2] that the condition (4.1) from Theorem 
8 (the growth of the variance) is satisfied for the Airy kernel (K t = K from 
(2.36), It expanding), the Bessel kernel (K t = K from (2.37), It expanding) 
and for the families of kernels {K n } corresponding to random matrices from 
the Classical Compact Groups (§2. 3b), §2. 3c)). In all these cases Var t #/ t 
growth logarithmically with respect to E t #/ t . 

Remark 13. To construct an example of the kernel < K < Id such that 
E#i- n ,n] = TrK • X[~n,n] ► oo as n ► oo, but Var#[_ Ttin] = Tr(K • X[-n,n] - 
(K ■ X[~n,n]) 2 ) stays bounded, consider {<p n {x)}n=-oc satisfying 

a) supp (p n E (n, n + 1), 

b) \\(Pn\\tf = 1- 

Then K(x,y) = E^-ool 1 ~ W+l) ' <Pn( x ) ' vM is the desired kernel. 

Indeed, E# [_ n , n] = ELJl - PTi) ™ °°' Var # = EL-nl 1 ~ PTi) ' 

PTT ^ ^"ool 1 " FTi) ' FTT < 00 • From the other side if < K < Id is 
compact, locally trace class and TiK ■ X[-n,n] +oo, then TrK ■ X[-n,n] — 
(K ■ X { _ nM f -> +oo. 
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The result of Theorem 8 can be generalized to a finite number of intervals. 

Namely, if . . . , 1+™^ are disjoint subsets such that Gov t (#jk), # T u))/Vt ~^ 

bij as t — > oo, 1 < i, j < m, where V t is some function of t growing to infinity, 

i 

then the distribution of ((# 7 (fe) — E t # J (fc) /V t 2 )k=i,..., m converges to the Tri- 
dimensional centralized normal vector with the covariance matrix (&ij)i<i,j<m 
(see [So2]). 

Finally, we turn our attention to the problem of the global distribution of 
spacings. Let E = W 1 or Z d , {Bj} k - =1 be some bounded measurable subsets 
of E, and {rij} k - =1 be some non- negative integers. We will be interested in 
the counting statistics of the following type 

rj L (B 1 ,...,B k ;n 1 ,...,n k ) := #(^ G [-L,L] d : # Xi + Bj =nj,j = l,...,k) 

(4.2) 

We can assume without loss of generality that {Bj} are disjoint and do 
not include the origin. If d — 1, k — 1, B x — (0, s], then r/ L ((0, s]), 0) is 
the number of the nearest spacings in [— L, L] greater than s : r/ L ((0, s], 0) = 
G [— L, L] : x i+ i — Xt > s}, and tjl((0, s]),n) is the number of n-spacings 
greater than s : t)l((0, s],n) = G [— L, L] : x i+n+ i — Xt > s}. In [Sol] 

we proved the convergence in law of the process UiA^fiAMz^SiJJ&fh^l ^ Q 

L2 

limiting Gaussian process in the case K(x,y) = s ^^. Recall that the 
convergence in law (Functional Central Limit Theorem) implies not only the 
convergence of the finite-dimensional distributions, but also the convergence 
of functionals continuous in the appropriate (e.g. locally uniform) topology 
on the space of sample paths. The proof of the Central Limit Theorem for 
the finite-dimensional distributions of ?7l((0, s], 0) can be extended essentially 
word by word to the case of arbitrary, not necessarily translation invariant, 
kernel K(x,y) and dimension d > 1, assuming the conditions (4.33), (4.34), 
(4.35) are satisfied. One can also replace (0, s] by an arbitrary measurable 
bounded B C E. For the convenience of the reader we sketch the main ideas 
of the proof of the finite-dimensional Central Limit Theorem below. Let 
us fix Bi, . . . , B k ; ni, . . . , n k . We construct a new (called modified) random 
point field such that i]l(Bi, . . . B k ; ni, . . . , n k ) is equal to the number of all 
particles of the modified random point field in [-L, L] d . Namely we keep 
only those particles of the original r.p.f. for which 

#x i +B j =n j ,j = l,...,k, (4.3) 
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and throw away the particles for which (4.3) is violated. The modified r.p.f. 
in general will no longer be a d. r.p.f. What is important is that its correlation 
functions and cluster functions (see Definition 6 below) can be expressed in 
terms of the correlation functions of the original determinantal r.p.f. Let us 
denote by pe(xi, . . . , xf, B±, . . . , B k ; n±, . . . , n k ) the £-point correlation func- 
tion of the modified r.p.f. Suppose that 

Xi&Xj + Bp, l<i^j<£, 1 < p < k. (4.4) 

Then by the inclusion-exclusion principle 

(-l) m 



mi 

m=0 



(4.5) 



ptfa, ...,x e ;B 1 ,...,B k ;ni,..., n k ) = ^ 

/ / 

< £ times > 

/ Pe+e-n+m(xi, . . . ,Xf, 

J((xi+uU B i) U- Ufo+U U B 5 ))m 

x u ,..., x ln , x 2 i, • • • , x 2n , xu, xen, yi, ■ ■ ■ y m )dyi . . . dy m 
dx a . . . dx in . . . dx u . . . dx ln , 
n = ni H h n k . 

If (4.4) is violated then the formula is quite similar, the only difference is 
that the exponent rij in (x* + Bj) n i = (x* + Bj) x • • • x (xi + Bj), 1 < % < 

< — rij times ► 

£, 1 < j < k, has to be replaced by rij — #(1 < r ^ i < k : x r E Xi + 
Bj). While formulas (4.5) appear to be cumbersome and lengthy, they are 
nevertheless quite useful for calculating the asymptotics of the moments of 
t)l(Bi, . . . , B k ; ni, . . . n k ). (Of course the assumption that the correlation 
functions of the original r.p.f. are the determinants is the key here.) Recall 
the definition of the cluster functions. 

Definition 6. The t-point cluster functions r^(xi, . . . , x^), £ — 1, 2, . . . , of a 
random point field are defined by the formula 

m 

r^, ...,x e ) = ^(-lr-V - 1)! • Up\gMGj)) (4-6) 

G j=l 

where the sum is over all partitions G of [£} = {1,2, ...,£} into subsets 
Gi, . . .,G m ,m = !,...,£, and x(Gj) = {x { : i e Gj}, \Gj\ = #(Gj). 
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The cluster functions are also known in the Statistical Mechanics as the 
truncated correlated funciton and the Ursell functions. Sometimes in the 
literature the r.h.s. of (4.6) defines (— l) l ~ 1 ri. Correlation functions can be 
obtained from cluster functions by the inversion formula 

m 

p e ( Xl , ...,x e ) = J2H r \ G M G j))- ( 4 - 7 ) 

G j=l 

((4.6) is just the Mobius inversion formula to (4.7).) The integrals of 
rg(xi, . . . , xe) over [-L, L] d x • • • x [— L, L] d = [-L, L] ld are closely related 

< 1 times > 

to the cumulants Cj(L) of the number of particles in [— L, L] d : 
V 1 (L) = J [ _ L>L]d r 1 (x 1 )dx 1 = Ci(L) = E# [ _ L)L]d , 

V 2 (L) := J [ _ LL]d J [ _ LL]d r 2 (x 1 ,x 2 )dx 1 

dx 2 = C 2 (L) - d(L) = Var# hLjL]d - E# { _ L<L]d , 

V s(L) := I[-L,L] d I[-L,L] d I[-L,L] d 

r 3 (xi, x 2 , x 3 )dxidx 2 dx 3 = C 3 (L) - 3C 2 (L) + 2C\(L). 
In general, 

n=l n=l 

(see [CL], [Sol]). For the determinantal random point fields 

rg(x 1 ,...,xe) = (-l)*' 1 ^2 K(x 1 ,x 2 ) • K(x 2 ,x 3 ) K(x t ,xi), 

cyclic cr£Si 

(4.9) 

where the sum in (4.9) is over all cyclic permutations, and the term written 
in the body of the sum corresponds to a — (1 2 3 . . .£). One can also rewrite 
(4.9) as 

r t {x 1 ,...,xt) = {-if' 1 ■ - B y2K(x aW ,x a{2) ) ■ 

t t^ l (4.10) 

K(x a ( 2 ),x a ( 3 )) K{x a ^),x a (i)). 

We note that the difference between (4.9) and the formula for £-point corre- 
lation 

p e (xi, ...,x e ) = ^(-1)^(0:1,^(1)) • K(x 2 ,x a{2) ) ■ K(x,x a(e) ) (4.11) 
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is that the summation in (4.9) is only over cyclic permutations. It appears 
that a relation between pg(xi, . . . ,xf, Bi, . . . , B k ; ri\, . . . , n k ) and rg(xi, . . . ,xf, 
Bi, . . . , Bk] ni, ... , nk) (at least when (3.19) is satisfied) is of a similar nature. 



Lemma 7. Let (4-4) be satisfied. Then 

-iy 



r e (xi, ...,xf,Bi,...,B k ;m,..., n k ) = ^ 



mi 

m=0 



-I times — ► (4-12) 



/ -I 

^(xi+Si^ix-xCxi+Bfc)"* J(x 
< — £ times — 

(xi, . . . , Xg, Xn , . . . , Xi n , 

■/((*i+U*=iBi)U-U(^+U*=i^)'» 
x 2 i, ■ ■ ■ , x 2n ■ ■ ■ , %, • • • , x £n , i/i, ... , y m )dyi . . . dy m 
dxgi . . . dxg n . . . dxn . . . dxi n , 

where pg + g. n+m / is defined below in (4-13). 

To define Pi+i. n +m,e recall that pg + g. n+m (xi, ...,y m ) = Y.aes l+l . n+m (~ l Y 
K{xi,a{xi)) • ... ■ K(y m ,a(y m )), where a is a permutation on the set of 
variables (x ± , . . . , x e , x u , . . . , x ln , y x ,...y rn ). We write 

(-\yK(x 1 , a(xi)) ■ K(y m , a(y m )), 

•n+m 

(4.13) 

where the summation in i s over the permutations a satisfying the follow- 
ing property: 

Let r be a multivalued map defined on {1,...,£} with the values in 
{!,...,£}: 



r(i) = {j : a^{xi,x iU . . . ,x in }\^^{y u . . . ,y m }^\(xi + [J B p ))^j 



(4.14) 
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then for any 1 < i, j < £ there exists N = N(i,j) such that 

r N {i)3j. (4.15) 

Remark 14. The proof of Lemma 7 in the case d = l,K(x,y) = ^I'f , 
B\ = (0, s],rii = 0, was given in §3 of [Sol]. In the general case the argument 
is absolutely the same. As a corollary of Lemma 7 we obtain 

Lemma 8. Let 

\K(x,y)\<il>{x-y), (4.16) 

and (4-4) hold for the t-twple (x 1 , . . . , xg). Then for any 5 > the following 
estimate takes place: 



\re(xi, . . . , xe; . . . , B k ;n!, . . . ,n k \ < const(£, 5) ■ 
^ f ip(x 2 - xi) ■ ip(x 3 - x 2 ) • ■ ■ ■ ■ ip(xi - x e ) J 



(4.17) 



cyclic <r£Se 



□ For the proof of Lemma 8 we refer the reader to [Sol] §3. The key element 
of the proof is the upper bound on the absolute value of the m th term in 
(4.12) by 

consti(n, £) ■ — -const™ • mini const3(n, £);(£ + in + m)\ ■ 
ml I 

^ f *P(x2 - xx) ■ ip(x 3 - x 2 ) ■ ■ ■ ■ ■ ip(xi -xe)\ >. 



cyclic a&S( 



□ 

1-5 



If ip 1 6 G L 2 (E) for some < 5 < 1, then Jj_ L L ^ d - ■ ■ Jj_ L L ^ d ip(x 2 — X\) 
. . . • ip(xi — xg) 1 ~ s dxi . . . dxe < const(^) • Jj_ L £ j d i>(x — y) 2 ~ 2S dxdy = 0(L d ), 
therefore by Lemma 8 

/ r i (x 1 ,...,xf,B 1 ,...,B k ;n 1 ,...n k ) 

J[-L,L] M n(4A) (4.18) 

dx!...dx t = 0(L d ),£= 1,2,... 
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In particular 



Erj L (B 1 , . . . , B k ; n u . . . , n k ) = Vy(L) = 
/ ri(x; B 1} . . . ,.B fe ;ni, . . .,n k )dx = 0(L d ) 

J[-L,L] d 

Suppose that one could show 

Var 77 L (Si, . . . , B k ; m, . . . n k ) = V X {L) + V 2 (L) = 
/ r 1 (x;B 1 ,...B k ;n 1 ,...,n k )dx+ / / 

J[-L,L] d J[-L,L] d J[-L, 

r 2 (xi, x 2 ; By,..., B k ; m, . . . n k )dxidx 2 = const • L d (l + o(l)), 



(4.19) 



(4.20) 

L] d 



/ r e (x u ...,x e ;B 1 ,...,B k ;n 1 ,... n k )dxi ...dx t 

J[-L,LY d \(AA) ( 4 21) 

= o(Lf),£>2. 

Since the I th cumulant of t]l is a linear combination of Vi(L),i = 1,2, . . .1 
(see (4.8)), the estimates (4.18)-(4.21) would imly that the £ th cumulant of 

Id 

t]l is const • L ■ (1 + o(l)) for £ = 2 and grows slower than L 2 for I > 2. 
This in turn would imply that while the second cumulant of ^==^= is 1, all 
the other cumulants of ^~ Er?I - go to zero as L — > +00. The last statement 

is equivalent to the statement that the moments of ^gft= converge to the 
moments of the normal distribution, and in particular 

3^-=.JNr<o,i). 

VVar r] L 

Of course the devil is in the details. It turns out that there is no nice extension 
of the formulas (4.12), (4.13) to the case when (4.4) is not satisfied. Below 
we show how one can overcome these difficulties in the case of rj L (B; 0) (i.e. 
k — 1, rii = 0). We introduce the centralized £-point correlation functions by 
the formula 



m 



p ( f\x 1 ,...,x e )=Y^ IIW^Gi), ( 4 - 22 ) 

G j=l 
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where Y^** is the sum over all partitions G = {G±, . . . , G m }, m = 1, 2, . . . of 
{1, . . . ,£} into two- and more element subsets (i.e. \Gj\ > l,j = 1, . . . ,m). 
It follows from (4.7), (4.22) that 



\x u ...,x t ) = p e (x u ...,x e ) + ^2(-l) p II 

p=l l<h<—<i p <£ s=l 



Pi(xi.) • pt- P ((xi,...,xe) \ (x h ,. . . ,x ip )^j = p t (x u ...,xt) - 

t P s \ 

Yl *[Ipi( x u) 'Pe-p[ {xi,...,x t ) \ (x h ,...,x ip ) J. 
p=i i<h<-<i p <e s=i ^ ' 



(4.23) 



Let us denote by M^(L) the integral of the centralized £-point correlation 
function of the modified random point field over [-L, L] M , 

M$(L) = [ ■ ■ ■ [ pf(x u . . . , xt- B x - 0)d Xl . . . dx e . (4.24) 

J{-L,L] d J[~L,L] d 

We have 

• - 1) • . . . • (r) L - i + 1) = e- m ^ • e^" 1 )^ • (4.25) 

If we can show that 

, , f (2n - 1)!! • const? • L nd • (1 + o(l)) for £ = 2n, 



and 

E77 L = const 2 -L d -(1 + 0(1)), (4.27) 

then (4.25) implies 

. f (2ra - 1)!! • (const! + const 2 ) n • L nd • (1 + o(l)) for £ = 2n, 
E ^- E ^ = {o (l¥) f 0r « = 2„ + l, 

(4.28) 
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and 

7] L - Er] L w 

2 > N(0, consti + const 2 ) 

One can in principle calculate M^ C \L) from (4.12), (4.13). Indeed, if 

Xi - xj B, (4.29) 

(we remark that (4.29) is exactly (4.4) written in the case k — 1, n± — 0), then 
the expression for pf\x\, . . . , xf, B; 0) can be obtained from (4.22), (4.12), 
(4.13). Otherwise pe(xi, . . . ,xf, B;0) = 0, and (4.23) implies 

1 V 

P { ;\ Xl ,...,xf,B;0) = J2(-lT- E II 

p=i i<h<-<i p <e s=i (4.30) 

ri(x is ;B;0) ■ p e _ p ([x u . . . , x e ) \ (x h , . . . , x ip )^j . 

If for an (£ — p)-tuple (x±, . . . , xt) \ (x^, . . . , Xi p ) the condition (4.29) is not 
satisfied, then the corresponding term pe_ p ((xi, . . . , xg) \ (x^, . . . , Xi p )) in 
(4.20) is zero. If (4.29) is satisfied for (x± : . . . ,xe) \ (x ii: . . . ,x ip ), then we 
iterate (4.23) again 

Pe-pl(x 1 ,...,x e ) \ {x h ,...,x ip ) 



We claim 



Lemma 9. Let the condition (4-29) be not satisfied for the £-tuple (x 1 , . . . , X(). 
Then 



x 1 ,...,x t ;B;0)= ^ ^ • JJn^; S; 0) • (4.31) 

$CDC{1,...,£} 



where 



c D = £ (- 1 ) |A| - ( 4 - 32 ) 

ADD, 

(4-29) is satisfied for x(A) 
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In particular Co = if (4-29) is not satisfied for x(D) or if there exists 
1 < % < £, % ^ D, such that for any 1 < j < £ Xi — x 3 ■ (jL B U (-B) . 
Proof easily follows from the above arguments. 

Theorem 9. Let (X, B, P) be a determinantal random point field with the 
kernel 

\K(x,y)\<^(x-y), (4.33) 

where ip is a bounded non-negative function such that ip- (log(^^)) n G L 2 (E) 
for any n > 0. Let for i] L (B; 0) = G [-L, L] d : + B) = 0) we have 

Var ^(5; 0) = a 2 ■ L d -(1 + 0(1)) (4.34) 

TTien t/ie Central Limit Theorem holds: 

ri L (B,0)-Eri L (B;0) ^ ^ 
L? 

Remark 15. If Cov^^; 0), ^{Bf, 0)) = % • L d • (1 + 0(1)), l<i,j<p, 
then 

^(gt;0)-E^(^;0) \ „ 

-5 ^(0>(o«)l<ij<p)- I 4 - 35 ) 

L2 J l<i<p 

Recall that 

CovMBf, 0); n L {Bj- 0)) = E( VL (B t ; 0) - Et^, 0)) • (rj^Bf, 0) - 

^□{xi-xa^BiUC-Sj)} m=0 m - ^((xi+B i )U(x2+B J -))"> 

P2+m,2(a:i, ^2; yi, • • • , ym)^yi • • • dy m )dx l dx 2 ~ / ri(xi; 5*; 0) • 

J[-L,L] d 

/ ri(x 2 ; Bj-,0)dx2dxi + / ri(x; £>i U -B2, 0)dx 

V(xi+B i )U(xi-B J -) J[~L,L] d 

Remark 16. Lemma 8 suggests slightly more restrictive condition on ■?/>, 
namely vp 1 ' 6 G L 2 (E) for some < 5 < 1. However, looking at the proof 
of Lemma 6 one immediately realizes that it is possible to replace ip 1-6 in 
(4.17) by i) ■ (log(^i))™ with n > U. 
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Proof of Theorem 9. It follows from (4.25)-(4.28) that it is enough to 
show 

/ p£(* ^o)^...^ 

J[-L,L] 2 " d 

(2n-l)!!-(/ { _ LM M n r 2 (x,y;B;0)dxdy - 

\J {x-ygBn(-B)} 

;B;0) f n(y;B;0)dydx) +d(L nd ), 

J(x+B)U(x-B) J 



\-l,lY (4.36) 



T\\X 

1 (x+B)U{x-B) 
71 = 1,2,..., 



/(c) \ - ( 2n + l A 

Pk+iv^ 1 ' • • • ' x 2n+i; -B; 0)dxi . . . dx 2n +\ — o[L 2 
_ LiL] (2n+l)d V / (4.37) 

n = l,2,... 



□ 



Lemma 10. 



/ Pzlixu ■■■■> x 2n, B; 0)dxi . . . dx 2n = 

J[-L,L] 2nd n(4.29) 

(2n-l)\\-([ { _ L ^ dn r 2 (x,y;B;0)dxdy) +d(L nd ), 

\J{x-ygBU(-B)} J 



(4.38) 



/(c) \ —/ 2n + ! (t\ 

p 2r i +1 {x 1 ,...,x 2n+1 ;B;0)dx 1 ...dx 2n+1 = o(L 2 ). 
1 -L,L]( 2 "+ 1 ) d n(4.29) 

(4.39) 

Recall that all r e (xi, . . . , xg; B; 0) are bounded functions (see (4.17)). Let 
us rewrite (4.22) as 

m m 
G j=l G j=l 

where Yl' i s the sum over all partitions of {1, . . . ,£} into pairs, and Y" is 
the sum over all other two- and more element partitions. Let £ be even, 
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i = 2n. Integrating YI'g over [—L,L] 2nd fl (4.29) we obtain exactly the 
r.h.s. of (4.38) (there are (2n — 1)!! partitions of {l,...,2n} into two- 
element sets). It follows from (4.17) and the estimate below Lemma 8 that 
f, L L ud \rt{xi, . . . , xe; B; 0)\dxi . . . dxe = 0(L d ). Therefore the integral of 

Y2g over [-L,L] 2nd n (4.29) is d(L nd ). The formula (4.39) can be proven in 
the same way. □ 
To estimate 



p ( 2n{xi, ■ ■ -,x 2n ; B; 0)dxi . ..dx 2n (4.40) 

L,L]2nd\( 4 29) 



we introduce the equivalence relation on {x±, . . . , x 2n } by calling 
bors" if there exists a sequence of indices 1 < iq, ii, . . . , i u < 2n, 1 < u < 2n, 
such that io — i, i u — j, and Xi s+1 — Xi a G B U (-B), s — 0, . . . , u — 1. We 
claim that the contributions of order 0(L nd ) appear in (4.40) only from such 
sets of (xi, . . . x 2n ) where each equivalence class of "neighbors" has either one 
or two indices. Consider for example the case when we have k two-element 
classes {xi,x 2 }, . . . , {x 2 fc_i, x 2 k} and 2n — 2k one-element equivalence classes 
{x 2 k+i}, • • • , {x 2n }- Similarly to the calculations on pp. 596-597 of [Sol] we 
verify that the integral of /4n (x\, . . . , x 2n ; B; 0) over the subset of [— L, L] 2nd 
corresponding to the above partition is equal to 

(2n-2k-l)\\- (- [ n^SjO) / n(y; B;0)dydx 

V J{-L,L] d J(x+B)U(x-B) 



(/ 



-k 



r 2 (x,y;B;0)dxdy) +d{L nd ). 
\-L,L] 2n r\ J 

{x-^flnl-B)) 

(4.41) 

After the summation over all partitions into one- and two-element equivalence 
classes of "neighbors" (we remark that (4.38) corresponds to the partition 
into singletons), we obtain exactly (4.36). It follows from Lemma 7 and 
(4.17) that all other partitions into the equivalence classes give negligible 
contributions. (4.37) can be proven in a similar fashion. □ 

The conditions of Theorem 9 are very unrestrictive in the case of trans- 
lation invariant kernels. The covariance function of the limiting Gaussian 
process w. lim ^((°' s ]' H^((M>o) j g ^ en gj ven by the (i-dimensional ana- 
logues of the formulas (37), (38), (26) from [Sol] (of course one has to re- 
place S1 °^~^ by K(x — y)). Here and below we denote by (0, s] the rectangle 
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(0, Si] x • • • x (0, Sd], s — (si, . . . , Sd). In particular, if K{x) is continuously 
differentiable the limiting Gaussian process is Holder-continuous with any 
exponent less than |. Among other characteristics of the modified random 
point field (with respect to B = (0, s], n = 0) one may be interested in 
the spectral measure of the restriction of the group {£/*} to the subspace 
of the centralized linear statistics. We shall denote the spectral measure by 
fj>( s \d\). Recall that the spectral measure fx^(dX) = n(dX) of the original 
determinantal random point field is given by (3.12). In particular, for the 
sine kernel 

dfj, 
dX 

After lengthy, but rather straightforward calculations one can obtain that in 
the case of the sine kernel : 

<*//•> *V , |A| / 4,3, , , „, m , „ m2l , 




dX 9 +£-|i-3-V|+o<**)+o(WV)+o(|V|.o 



(4.42) 



We note that -^-(0) ^ if s ^ 0, s small, which is consistent with Var 
t)l({0, s];0) ~ L. For the proof of the Functional Central Limit Theorem we 
refer the reader to pp. 577, 598-600 of [Sol]. Suppose that 

L- d —r] L ((0, s]; 0)), L~ d — Cov(77 L ((0, s}; 0); ^((0, t\; 0)) (4.43) 

are uniformly bounded in L, s, t, where s, i belong to compact subsets of 
R^(Z+). By smoothing with a C°° approximate 5-function one can construct 
a continuous approximation fj L ((0, s]; 0) such that |?7l((0, s]; 0)— r/ L ((0, s] ; 0) j < 
1. As a result 

7y L ((0,s];0)-E7y L ((0,s];0) 
r ^ 

L2 



is a random continuous function in s, and 

77 L ((0, s]; 0) - Ef) L ((0, s]; 0) rj L ((0 : s]; 0) - E^((0, s]; 0) 



r- r- 

L2 L2 



<4 (4-44) 

L2 



The distribution of the random process ^((M;0)-e^((o,s];0) defines a prob _ 

ability measure on C([0, oo) d ). By the convergence in law of random pro- 
cesses we mean the weak convergence of the induced probability measures 
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on C([0, oo) d ) (see [B], in general one can consider different spaces of sample 
paths, e.g. the space of cadlag functions, instead of the space of continuous 
functions). 

Theorem 10. Let the condition (4.33), (4-34), (4-35), (443) be satisfied. 
Then the random process 

iy £ ((0,s];0)-E^((0,fl;0) 
Lh 

converges in law to the limiting Gaussian process. 
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